0 Oh  ‘■]o‘y 


AD-A247  554 


PL-TR-91-2206 


INVESTIGATION  OF  RING  CURRENT  /  STORM  DYNAMICS 


H.  L.  Collin 
J.  B.  Cladis 
J.  M.  Quinn 


Lockheed  Palo  Alto  Research  Laboratory 
Space  Sciences  Laboratory 
3251  Hanover  Street 
Palo  Alto,  CA  94304-1191 


August  31  1991 


./f 


a;/ 


biacl'. 


aud 


Scientific  Report  No.  1 

Approved  for  public  release;  distribution  unlimited 


PHILLIPS  LABORATORY 

AIR  FORCE  SYSTEMS  COMMAND 

HAN3COM  AIR  FORCE  BASE,  MASSACHUSETTS  01731-5000 


92-04155 


This  technical  report  has  been  reviewed  and  is  approved  for  publication 


MICHAEL  D.  VIOLET,  LT,  USAF 
Contract  Manager 


E.  G.  MULLEN 
Branch  Chief 


^  (■< 


Y.' 


WILLIAM  SWIDER 
Deputy  Director 


This  document  has  been  reviewed  by  the  ESD  Public  Affairs  Office  (PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS). 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical 
Information  Center.  All  others  should  apply  to  the  National  Technical 
Information  Service. 


If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing  list, 
or  if  the  addressee  is  no  longer  employed  by  your  organization,  please  notify 
GL/I'IA,  Hanscom  AFB,  MA  01731.  This  will  assist  us  in  maintaining  a  current 
mailing  list. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices  on 
a  specific  document  requires  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No  0704  0188 


Public  ceporting  burden  (cc  ilii^  coiifclion  of  inforfriatiof  eiliffiaied  to  atfetoqe  '  hour  per  recporse  including  the  tune  tor  reviewing  instructions  searching  enstmg  data  sources, 
gathering  and  maintaining  the  data  needed  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  an,  other  aspect  of  this 
collection  of  information,  including  suggestions  lor  routing  this  burden  to  UVashington  Headquarters  Services  Directorate  tor  information  Operations  and  Reports,  i  ^  t  S  Jefferson 
Oavis Highway. Suite  1204  Arlington  UA  22202-4302.  and  to  the  Office  of  Management  and  Budget  Paperwork  Reduction  Proiect(0704-0l88).  Washington.  DC  20503 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  OATES  COVERED 


8/31/91 


4.  TITLE  AND  SUBTITLE 


Investigation  of  Ring  Current/Storm  Dynamics 


6.  AUTHOR(S) 

H.L.  Collin,  J.B.  Cladis,  J.M.  Quinn 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES) 

Lockheed  Palo  Alto  Research  Laboratory 
Space  Sciences  Laboratory 
3251  Hanover  Street 
Palo  Alto,  CA  94304-1191 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOOR£SS(ES) 

Phillips  Laboratory 
Hanscom  APB,  MA  01731-5000 

Contract  Manager:  Michael  Violet  Lt,  USAF/PHP 


S.  FUNDING  NUMBERS 

P19628-90-C-0119 

PE62101F 

PR7601 

TA22 

WVLC 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 

PL-TR-91-2206 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  unlimited 


13.  ABSTRACT  (Maximum  200  words) 


12b.  DISTH't;UTION  CODE 


This  report  outlines  the  achievements  of  the  first  year  of  work  on  data  acquired 
by  the  Lockheed  Ion  Mass  Spectrometer  (IMS-LO),  one  of  the  space  radiation  group 
of  instruments  on  the  Combined  Release  and  Radiation  Effects  Satellite  (CRRES). 
The  purpose  of  IMS-LO  is  to  study  the  low  energy  (llOeV-35  keV)  ion  component  of 
the  ring  current  plasma.  The  primary  data  processing  is  proceeding  rapidly,  as 
is  the  construction  of  databases  of  validated  low  energy  ion  data.  A  number  of 
events  have  been  selected  for  special  study  and  one  of  these  has  been  analysed 
in  detail.  A  tentative  identification  has  been  made  of  a  new  mechanism  which 
causes  loss  of  ions. 


14.  SUBJECT  TERMS 


Geomagnetic  Storms,  Ring  Current,  Ions,  Radiation  Belts 


15  NUMBER  OF  PAGES 

86 


16  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18  SECURITY  CLASSIFICATION  19  SECURITY  CLASSIFICATION  20  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 


Unclassified 


NSN  7540-0  I  ^SO-SSOO 


Unclassified 


Unclassified 


SAR 


Standard  form  298  (Rev  2  89) 


CONTENTS 


page 


1.  Introduction  1 

2.  Experiment  Status  and  Operations  3 

3.  Progress  and  Current  Activities  5 

3.1  Generation  of  Lockheed  IMS  data  base  5 

3.2  Empirical  model  of  ring  current  ion  composition  9 

3.3  Magnetospheric  dynamics  inferred  from  ion  transport  9 

3.3.1  Introduction  13 

3.3.2  Dynamics  of  the  Storm  of  26  August  1990  14 

3.3.2. 1  Observations  14 

3.3.2.2  Ion  Pressures  23 

3.3.2.3  Interpretation  of  Data  33 

3.3.2.3.1  Ion  Sources  33 

3.3.2.3.2  Variations  of  Particle  Presiures  34 

3.3.2.3.3  A  New  Loss  Mechanism  35 

4.  Future  plans  38 

References  39 

Appendix  1:  The  Low  Energy  Ion  Mass  Spectr  'ter  (ONR  307-8-1, fi)  41 

Appendix  2;  Abstrjicts  of  Papers  Presented  52 

Appendix  3:  Distribution  in  Magnetotail  of  Ions  56 


iii 


Acoasf 

;ton  For  j 

’"mtis 

DTIC 
Uiia’.i  ■ 

Ju 

GRA&I 

CV3  □ 

■'.iiiced  0  j 

uy - - - -  - - 

m  stril-'Jtion/ 

Avs  • 

lability  Codes 

Dlat 

nil 

jAVail  and/or 

1  Spacial 

1.  Introduction 


The  effort  consists  of  three  Tasks  in  support  of  one  of  the  major  goals  of  the  Combined 
Release  and  Radiation  Effects  Satellite  (CRRES)  which  is  to  develop  a  time-dependent 
model  of  the  magnetospheric  radiation  environment  that  presents  a  hazard  to  the  perfor¬ 
mance  of  spacecraft  electronic  systems.  This  effort  will  model  the  low  energy  ion  species, 
H+,  0+,  He+,  and  He++. 

There  is  no  question  regarding  the  need  for  such  modeling.  These  heavy-ion  species  - 
especially  the  ones  in  the  ring  current  with  energies  extending  to  300  keV/e  range  [Gloeckler 
et  al.,  1985]  -  contribute  appreciably  to  the  degredation  of  spacecraft  thermal  control 
surfaces,  optical  coatings  and  filters,  and  other  sensitive  surfaces.  Moreover,  these  ions 
are  fundamental  to  the  dynamic  modeling  of  the  radiation  belt:  (1)  the  low  energy  ions 
strongly  effect  the  dynamics  of  the  magnetosphere  at  dist2inces  greater  than  about  7  Re 
in  the  magnetotail,  and  they  become  strongly  affected  by  the  magnetospheric  dynamics 
inward  of  about  6  Re  -  Thereby  offering  a  means  of  modeling  the  stormtime  geoelectric 
and  geomagnetic  fields  and  their  effects  on  all  the  trapped  radiation;  (2)  these  ions  are  the 
principal  source  which  is  eventually  accelerated  to  become  much  higher  energy  particles 
in  the  inner  radiation  belt;  and  (3)  these  ions  strongly  affect  the  redistribution  of  inner 
belt  ions  -  in  species  type  as  well  as  spacially  and  temporally  -  and  loss  through  mviltiple 
charge-exchange  processes  and  multi-species  instabilities  and  wave-paxticle  interactions. 

For  example,  Wrenn  (1989]  has  recently  reported  that  the  rate  of  recovery  of  magnetic 
storms  is  related  to  solar  cycle.  He  attributed  the  change  in  recoveiy  rate  to  changes  in  the 
ion  composition  which  strongly  infiuences  the  rate  of  charge  exchange,  one  of  the  principal 
loss  processes.  Another  important  ring  current  loss  process  is  pitch  angle  diffusion  by  ion- 
cyclotron  waves.  Kozyra  et  al.  [1984]  showed  that  the  generation  of  these  waves  by  the 
interaction  of  plasmaspheric  plasmas  with  the  ring  current  is  strongly  dependent  on  the 
composition,  densities  and  temperatures  of  the  interacting  plasmas. 

The  development  of  a  comprehensive  and  accurate  model  of  the  nng  current  and 
radiation  belts  requires  a  thorough  understanding  of  ring  current  dynamics,  which  in  turn 
reqmres  detailed  knowledge  of  the  particle  source  contnbutions,  transport  phenomena, 
energization  processes,  and  scattering  and  loss  mechanisms. 

The  low  energy  plasma  plays  a  major  role  in  the  dynamics  of  the  inner  magnetosphere. 
A  critical  component  of  this  plasma  is  the  composition  of  the  ion  component.  These 
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low  energy  ions,  which  originate  in  the  solar  wind  and  the  ionosphere  serve  as  a  source 
for  the  ring  current  and  radiation  belts.  The  principle  tracers  of  plasma  origin  in  the 
magnetosphere  are  measurements  of  ion  flux  ratios  and  charge  state  density  ratios. 

Because  of  limitations  in  available  satellite  instrumentation  measurement  of  ion  com¬ 
position  in  the  ring  current  and  radiation  belts  has  begun  only  fairly  recently.  The  com¬ 
prehensive  data  required  for  an  environmental  modal  were  not  available  before  CRRES 
[Williams,  1987]  and  earlier  efforts  at  modeling  the  behavior  of  the  low  energy  ions  were 
imperfect,  as  evidenced  by  the  number  of  unanswered  question  that  they  raised  [Delcourt 
et  al.,  1989).  A  principal  objective  of  this  experiment  is  to  provide  a  data  base  suitable  for 
modeling  efforts. 

Task  1  is  to  process  the  raw  Lockheed  ONR-307-8  Ion  Mass  Spectrometer  (IMS)  data 
from  the  agency  tapes  to  provide  the  basic  data  base  of  quality  controlled  calibrated  in  situ 
measurements  of  ion  composition  upon  which  the  scientific  studies  are  based.  The  other  two 
Tasks  make  use  of  this  data  base.  Task  2  is  to  construct  an  empirical,  statistical  model  of 
the  plasma  composition  and  distribution  function  as  a  function  of  magnetospheric  location 
for  quiet  and  disturbed  magnetospheric  conditions.  During  magnetic  storms  the  model  is  to 
be  organized  by  storm  phase  in  order  to  reflect  the  development  of  the  ring  current  plasma. 
Task  3  is  to  investigate  aspects  of  magnetospheric  dynamics  by  making  a  specialized  study 
of  magnetic  storm  events  selected  from  the  Lockheed  IMS  data  base.  For  each  event  it  will: 
1)  Determine  time  dependent  models  of  geomagnetic  and  geoelectric  fields  consistent  with 
the  transport  of  the  measured  ions  between  observation  points  on  successive  legs  of  the 
CRRES  orbit.  2)  The  transport  of  ions  from  typical  sources  of  up-flowing  ionospheric  ions 
will  be  computed  and  compared  with  Wk  CRRES  in  situ  measurements  of  magnetospheric 
ions.  The  ion  trsinsport  will  eilso  be  compared  with  possible  storm  predictors. 

We  will  make  use  of  observations  of  ion  species,  principally  ,  0"^,  He'^'^  and  He'*',  in 
the  energy  range  100  eV/e  to  about  100  keV/e.  These  ions  become  strongly  affected  by  the 
magnetospheric  dynamics  thereby  offering  a  means  of  modeling  the  stormtime  geoelectric 
and  geomagnetic  fields  and  their  effects  on  all  the  trapped  radiation.  They  are  the  principal 
soiuce  of  ions  which  are  eventueilly  accelerated  to  become  much  higher  energy  particles 
in  the  inner  radiation  belt  and  they  strongly  affect  the  redistribution  and  loss  of  inner 
belt  ions  through  multiple  charge-exchange  processes  and  multi-species  instabilities  and 
wave-particle  interactions. 

To  achieve  these  objectives  ONR  307-8-1,2  (IMS-LO)  is  designed  to  measure  plasmas 
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that  are  the  sources  of  radiation  belt  particles,  and  to  provide  data  on  the  origin  and 
acceleration  of  these  plasmas.  To  achieve  these  objectives  IMS-LO  measures  energy  and 
mass  spectra  covering  the  ranges  of  £/q  from  O.I  to  35  keV/e  and  M/q  from  1  to  32 
AMU/e  with  good  coverage  of  pitch  angles  throughout  the  orbit.  A  full  description  of  this 
instrument  is  in  Appendix  1 


2.  Experiment  Status  and  Operations 

CRRES  was  successfully  launched  on  25  July  1990  and  initialization  operations  of  the 
ONR  307-8  instruments,  which  provide  data  that  are  used  for  this  contract,  began  shortly 
thereafter.  Notable  events  in  the  initialization  sequence  are  sununarized  in  table  1. 


Table  1 

ONR  307-8-1,2  initialization  events 

Date 

Time  (UT) 

Orbit 

Event 

26  July 

21:40 

8.1 

Power  on 

29  July 

07:40 

9.1 

Memory  load  and  check 

2  August 

10:50 

19.1 

8-1  high  voltage  on 

2  August 

21:40 

20.1 

8-2  high  voltage  on 

4  August 

12:00 

24.1 

Configure  for  normal  operations 

26  August 

17:10 

78.1 

Reconfigure  for  SWEEP-LOCK  mode 

Initialization  preceded  smoothly.  Following  initialization,  ONR  307-8-1  eind  ONR  307- 
8-2  were  operated  in  SWEEP/LOCK  mode  to  give  good  sensitivity  and  temporal  2md  pitch- 
angle  resolution  of  -.e  major  species  H'*',  and  0"^,  while  retaining  full  mass  resolution. 

Both  instruments  continue  to  work  perfectly.  Almost  the  only  non-scheduled  loss  of 
data  occurs  as  a  result  of  the  loss  of  synchronization  caused  by  the  dropouts  in  the  vehicle 
clock.  When  these  occur  the  instruments  lose  synchronization  but  are  able  to  recover 
within  a  couple  of  minutes. 

The  only  other  anomaly  which  has  occured  to  date,  happened  on  18  July  1991  when 
the  CSTC  controler  reported  that  ONR  307-8-1  was  cooling  rapidly  and  its  high  voltage 
monitors  were  showing  abnormally  low  values.  The  status  of  the  instrument  was  examined 


3 


on  the  next  available  contact  and  it  was  ascertained  that  one  of  its  control  registers  had 
become  incorrectly  set  and  that  in  consequence  two  of  the  high  voltage  generators  had 
been  di6>enabled.  The  fall  of  instrument  temperature  was  the  result  of  the  reduction  in 
power  consumption.  Commands  were  sent  at  once  to  reset  the  control  register  and  restart 
the  high  voltage  generators.  The  precise  cause  of  the  registers  becoxmng  unset  is  not  clear. 
Geophysical  conditions  were  very  disturbed  at  the  time  and  it  is  probable  that  a  charging 
event  or  SEU  was  the  ultimate  cause. 

The  failure  of  battery  2,  on  22  December  1990  led  to  duty  cycling  of  the  ONR  307*8 
instruments  during  eclipses  in  order  to  conserve  the  limited  battery  power.  The  instrument 
temperatures  were  monitored  and  although  they  cooled  considerably  they  did  not  reach 
dangerously  low  levels. 

Data  analysis  has  been  proceding  steadily  and  all  the  scientists  ha”e  taken  part  in  CE¬ 
RES  science  working  groups.  The  American  Geophysical  Union  devoted  a  special  session 
to  preliminary  reports  of  CRRES  data  at  its  Spring  Meeting  in  Baltimore.  Two  papers 
were  presented;  ‘Low  Energy  Ring  Current  Ions  during  the  Storm  of  26  August  1990;  Ob¬ 
servations  by  IMS-LO  on  CERES’,  by  H.  L.  Collin,  J.  M.  Quinn  J.  B.  Cladis,  J.  B.  Cladis, 
and  ‘Storm- Time  Dynamics  Implied  by  the  CRRES  Measurements  of  Ring- Current  Ions’, 
by  J.  B.  Cladis,  H.  L.  Collin,  and  W.  E.  Francis.  Abstracts  of  these  papers  are  in  Appendix 
2. 

In  addition  to  radiation  studies  an  objective  oi‘  the  CRRES  mission  is  a  series  of 
chemical  release  experiments.  The  releases  were  made  in  three  very  successful  campaigns. 
At  low  altitude  during  September  1990,  and  July  and  August  1991,  and  at  high  altitude 
during  January  and  February  1991.  ONR  307-8-1,2  supported  these  chemical  relras'e  cam¬ 
paigns  under  contracts  from  NASA.  During  the  high  altitude  chemical  release  campaign 
ONR  307-8-2  weis  run  in  operating  modes  specially  designed  to  enhance  its  sensitivity  to 
Ba'^  or  Li'*'  and  these  modes  were  commanded  in  real  time  to  match  the  release  selected  by 
the  chemical  release  Principal  Investigators.  Support  was  also  provided  by  supplying  the 
Principal  Investigators  with  real  time  data  on  the  plasma  conditions  in  the  magnetosphere. 
This  information  was  helpfxil  in  choosing  the  optimum  release  time.  During  the  low  alti¬ 
tude  champaigns  ONR  307-8-2  was  run  in  a  mode  sensitive  to  Ba"*".  Some  of  this  activity 
was  reported  in  a  paper  presented  at  the  Spring  Meeting  of  the  American  Geophysical 
Union  in  Baltimore  in  the  special  session  devoted  to  the  chemical  release  portion  of  the 
CRRES  mission;  ‘Real  Time  Monitoring  of  the  Geophysical  Psirameters  in  Support  of  the 
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CURES  Chemical  Releases’,  by  S.  B.  Mende,  E.  L.  Collin,  W.  K.  Peterson,  T.  A.  Fritz,  G. 
Baerendel,  R.  R.  Anderson,  D.  A.  Hardy,  H.  J.  Singer,  J.  C.  Foster,  and  D.  L.  Reasoner. 
An  abstract  of  this  paper  is  included  in  Appendix  2. 


3.  Progreis  and  Current  Activities 


3.1  Generation  of  Lockheed  IMS  data  base 

The  first  step  in  the  analysis  of  the  CRRES  data  involves  processing  the  agency  tapes 
provided  by  the  Phillips  Laboratory'.  The  agency  tapes  for  Lockheed  contain  data  from 
each  of  the  ONR  307  instruments  together  with  ephemeris  and  magnetometer  data.  A 
schematic  of  the  data  processing  sequence  showing  the  relationship  of  analysis,  databases 
and  data  products  is  in  Figure  1.  The  first  stage  of  proccessing  is  to  unpack  and  decode 
the  data  on  the  agency  tapes,  and  write  the  data  from  each  orbit  to  a  group  of  files 
corresponding  to  each  instrument  together  with  ephemeris  and  magnetometer  data.  At 
this  stage  all  data  is  in  the  form  of  counts  or  engineering  units  at  the  full  resolution  of  the 
instrument  and  is  independent  of  posri'wle  drifts  or  revisions  of  detector  calibrations.  This 
forms  the  bade  High  Resolution  Database  which  is  used  as  the  basis  of  all  other  analysis. 
It  is  held  on  a  series  of  optical  disks  each  of  which  holds  data  from  26  orbits. 

The  second  stage  of  data  processing  is  the  production  of  the  Survey  Plot  Library. 
This  consists  of  full  orbit  survey  plots  for  all  orbits.  This  library  forms  an  important 
intermediate  data  product  which  displays  IMS-LO  data  in  a  form  which  is  compact  and 
also  sufficiently  detailed  to  enable  periods  of  scientific  interest  to  be  selected  for  specialized 
study  under  Task  3  or  for  other  purposes.  The  survey  plots  are  also  an  important  tool 
for  the  assessment  of  instrument  performance.  For  all  survey  plots  the  data  are  averaged 
over  262s  intervals.  The  examples  shown  in  Figures  2  to  4  are  taken  from  the  storm  of 
26  August  1990  which  is  a  period  of  special  study  by  the  CRRES  Low  Energy  Pleisma 
Field  and  Wave  Working  Group  and  is  discussed  in  detail  in  section  3.3.2.  Figure  2  shows 
a  survey  plot  of  mass  spectra  as  a  function  of  time  from  orbit  77,  26  August  1990.  The 
upper  group  of  three  panels  show  mass  spectra  from  the  three  heads  of  IMS-LO-1.  Each 
panel  contaiins  horizontal  bands  near  the  bottom  and  top  which  correspond  to  and 
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0"^  respectively.  Fainter  discontinuous  bands  between  them  correspond  to  He"*"  and  He''"'’ . 
The  lower  group  of  three  panels  show  similar  mass  spectra  from  IMS-LO-2.  Survey  plots  of 
energy  spectra  and  pitch-angle  distributions  for  orbit  76  and  orbit  77  are  shown  in  Figure  3 
and  Figure  4  respectively.  The  upper  group  of  four  panels  show  energy  spectra  of  H''' ,  O''’ , 
He''’,  and  He''”'’.  For  these  plots  the  flux  has  been  sorted  into  15  logarithmically  spaced 
energy  bins.  The  lower  two  groups  of  three  panels  display  the  pitch-angle  distributions 
of  H'*’  and  0''  in  each  of  the  three  energy  bands  corresponding  to  the  three  heads  which 
make  up  the  IMS-LO  spectrometer.  For  each  of  these  the  flux  has  been  sorted  into  15* 
pitch-angle  bins. 

Special  event  studies  and  dynamic  modeling  which  need  maximum  resolution  and 
flexability  make  direct  use  of  the  High  Resolution  Databue.  Plots  of  energy  spectra  and 
pitch  angle  distributions  for  selected  periods  of  data  were  produced  as  part  of  the  study 
of  storm  of  26  August  1990  in  order  to  permit  more  detailed  examination  than  is  possible 
with  the  survey  plots.  Examples  are  shown  in  Figure  5  the  upper  panel  of  which  shows  O''’ 
energy  spectra  at  a  series  of  pitch  angles.  The  flux  equivalents  of  the  background  count 
rate  and  the  one  count  per  sample  are  shown  for  comparison.  In  the  lower  panel  the  same 
C’  data  displayed  as  pitch-angle  distributions  at  a  series  of  energies.  The  0''’  is  seen 
to  be  isotropic  at  energies  above  0.5  keV  while  at  lower  enerpes  it  displays  a  pitch-angle 
distribution  which  is  strongly  asymetric  with  a  marked  field  aligned  flow  at  pitch  angles 
above  120*.  These  distributions  are  discussed  further  in  section  3.3.2. 

Other  routine  uses  of  the  High  Resolution  Database  are  the  production  of  the  Sum¬ 
mary  Database,  which  is  described  in  section  3.2,  and  the  monitoring  of  instrument  status 
and  calibration.  Reviews  of  in  flight  calibration  has  ted  to  minor  revisions  of  some  of  the 
mass  settings  used  in  the  LOCK  mode.  These  revisions  resulted  in  a  modest,  but  worth¬ 
while,  increase  in  sensitivity.  We  have  made  a  provisional  alignment  of  the  calibration  of 
the  three  heads  which  make  up  each  of  the  ONR  307-8  instruments,  and  between  the  two 
instruments  ONR  307-8-1  and  ONR  307-8-2  which  enables  the  slight  differences  between 
their  geometry  factors  to  be  accounted  for. 

Initially  the  data  processing  was  carried  out  on  a  VAX-780.  Later,  however,  a  SUN 
system  was  made  available  and  the  major  part  of  the  data  processing  and  analysis  has 
been  transferred  to  this  system.  Eventually  almost  all  analysis  will  be  carried  out  on  this 
new  system.  The  transfer  of  software  caused  some  temporary  delay  which  was  rapidly 
overcome  due  to  the  much  greater  computing  power  of  the  SUN  system.  All  agency  tapes 
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rcdcved  Lo  date  have  now  been  processed. 


3.2  Empirical  model  of  ring-current  ion  composition 

The  preliminary  veraon  of  the  Empirical  Model  of  Low  Energy  Ring  Current  Ion 
Composition  will  consist  of  a  set  of  average  equatorial  energy  pitch-angle  distributions  of 
and  O'*" .  The  precession  of  the  CRRES  orbit  will  eventually  enable  the  model  to  cover 
all  magnetic  latitudes  as  well  as  L  values  up  to  about  7.  These  average  distributions  will 
be  compiled  from  ion  measurements  taken  along  the  CRRES  orbit  and  mapped  to  the 
magnetic  equator.  The  model  can  then  be  used  to  determine  the  average  ion  distributions 
away  from  the  equator  by  mapping  the  equatorial  distributions  down  the  field  lines  to 
the  location  of  interest.  As  data  accumulates  and  understanding  of  the  behaviour  of  the 
ring-current  plasma  improves  the  model  will  be  expanded  to  include  the  effects  of  magnetic 
storms. 

Increased  understzmding  of  the  behaviour  of  the  storm-time  magnetosphere  is  likely 
to  produce  improvements  in  the  models  used  to  describe  the  geomagnetic  field  and  so 
to  require  changes  in  the  mapping  of  ion  trajectories  between  CRRES  and  the  equator 
and  thus  revision  of  the  empirical  model  of  ring-current  ion  composition.  In  order  to 
avoid  the  computational  burden  of  having  to  generate  the  model  from  the  High  Resolution 
Database  each  time  a  revision  is  needed  an  intermediate  data  base  is  being  constructed. 
This  IMS-LO  Summary  Data  Base  is  a  compact,  partially  processed  data  base  consisting 
of  averages,  over  262  seconds,  of  IMS-LO  ion  and  electron  data  sorted  by  energy  and  pitch 
angle.  Data  from  both  IMS-LO-1  and  IMS-LO-2  is  combined.  Supporting  data  include 
time,  background  count  rate,  status  information,  measured  magnetic  field  and  ephemeris 
information.  This  Simimeiry  Data  Base  is  also  suitable  for  rapid  retrieval  of  energy  spectra, 
pitch  angle  distibutions  or  survey  plots  when  the  full  resolution  is  not  required,  and  as  a 
basis  for  statistical  studies.  The  ion  and  electron  data  are  recorded  as  count  rates  in  order 
to  madce  the  data  base  independent  of  revisions  to  the  instrument  calibration  which  cross 
calibration  studies  are  expected  to  indicate.  Calibrations  are  applied  to  the  count  rates 
when  the  data  base  is  accessed. 

3.3  Magnetospheric  dynamics  inferred  from  ion  transport 
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distributions  at  a  series  of  enersies. 


3.3.1  Introduction 


The  objectives  of  this  task  are  to  use  the  ion  data  during  disturbance  events  (substorms 
and  magnetic  storms)  to  infer  the  underlying  dynamics  of  the  magnetosphere,  e.g,  the 
time-dependent  magnetic  and  electric  fields,  and  sources  of  the  ions  injected  into  the 
trapping  region  of  the  magnetosphere.  The  ions  measured  with  the  IMS-LO  and  IMS- 
HI  instruments,  particularly  the  species  H'*’,  0"^,  He"*”^,  and  He'^  in  the  energy  range 
0.2  keV/e  to  about  300  keV /e,  will  be  used  in  this  study.  Further  information  useful  to  this 
study  include  the  onboard  magnetic  and  electric  field  measurements,  high  energy  proton 
and  electron  data,  IMP-8  solar-wind  measurements,  geomagnetic  indices,  and  available  ion 
data  from  other  spacecraft  (e.g.,  DMSP  and  Akebono). 

The  magnetic  and  electric  fields  will  be  studied  by  computing  the  drift  paths  (eastward 
and  westward)  of  the  ions  in  trial  held  models  between  one  leg  (inbound  or  outboimd)  of 
the  CURES  orbit  to  the  succeeding  leg.  Liouville’s  theorem  will  be  used  to  assess  the  trial 
field  models.  That  is,  if  the  fields  along  the  ion  paths  are  correct,  the  phase  space  densities 
of  the  ions  measured  at  the  satellite  at  the  two  points  will  be  the  same.  Corrections  wiU  be 
made  for  charge  exchange  losses.  Additional  ion-dependent  redistributions  or  losses  due 
to  instabilities  or  wave-particle  interactions  will  be  inferred  from  the  results  obtained  for 
the  different  species. 

One  of  the  complications  of  this  study  is  that  the  ions  are  not  only  affected  by  the  mag¬ 
netic  and  electric  fields,  but  they  also  affect  the  fields  through  their  collective  behaviour. 
During  substorms,  ions  of  energy  less  than  about  50  keV  appreciably  affect  the  magnetic 
and  electric  fields  in  the  near-earth  magnetotail  and  nightside  magnetosphere  (e.g.,  see  the 
attached  preprint  by  Cladis  and  Francis  [1991]  in  Appendix  3).  During  magnetic  storms, 
ions  of  energy  to  about  300  keV  affect  the  fields  throughout  the  magnetosphere.  Also, 
during  the  initial  phase  of  a  magnetic  storm,  the  fields  may  change  so  rapidly  in  time  and 
space  that  it  may  not  be  possible  to  obtain  reliable  estimates  of  the  vwation  of  the  fields 
along  the  ion  drift  paths.  For  these  reasons  the  study  will  begin  with  quiet-time  conditions 
and  proceed  toward  more  disturbed  times  as  experience  is  gained  on  the  evolution  of  the 
fields  and  the  sources  and  losses  of  the  ions. 

Good  data  on  the  low-energy  ion  fluxes  measured  with  the  IMS-LO  instruments  are 
already  available  for  the  study,  and  approximate  IMS-HI  fluxes  are  expected  to  be  available 
in  a  few  months.  The  ion  transport  models  and  preliminary  trial  magnetic  and  convection 
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electric  fields  are  already  available.  A  model  of  the  induced  electric  field  resulting  from 
compressions  and  rarefactions  of  the  magnetosphere  has  not  yet  been  developed. 

To  date,  the  principal  effort  on  Task  3  has  been  devoted  to  the  study  of  the  dynamics 
of  the  magnetic  storm  of  26  August  1990.  This  effort  is  described  below. 


3.3.2  Dynamics  of  the  Storm  of  26  August  1990 


3.3.2. 1  Observations 


The  IMS-8  measurements  of  the  interplanetary  magnetic  field  (IMF)  and  the  plasma  in 
the  solar  wind  were  kindly  sent  to  us  by  S.  Gussenhoven  of  PL/GD.  for  day  238  (26  August 
1990),  Leppings  measurements  of  the  total  magnetic  field  strength  and  its  components  are 
shown  in  Figures  6a  -  6d,  and  Lazarus’  measurements  of  the  plasma  density  and  speed  are 
shown  in  Figure  7.  The  significance  of  the  vertical  lines  in  Figure  7  will  be  discussed  later. 
Singer’s  magnetic  field  measurements  on  CRRES  in  VHD  coordinates,  which  he  designates 
as  X,Y,Z,  are  shown  in  Figures  8a  and  8b.  In  these  Figures  6Bx,6By,  and  SBg  are  the 
differences  of  the  measured  components  from  the  components  given  by  the  Olson-Pfitzer 
magnetic  field  model.  The  total  magnetic  field  strength  is  given  in  the  bottom  panel,  and 
the  broken-line  curve  in  this  panel  is  the  model  field  strength. 

Note  from  Figures  6  and  7  that  the  interplanetary  shock  is  first  observed  by  IMP- 
8  at  about  05:40UT  and,  from  Figure  8a,  that  the  sudden  commencement  is  observed 
on  CRRES  at  05:53UT.  At  the  time  of  the  sudden  commencement,  the  B^  component 
suddenly  increases  by  about  15  nT  and  long  trains  of  ULF  oscillations,  particularly  in 
the  By  component,  are  initiated.  At  about  06:15UT  the  interplanetary  Bg  component,  as 
shown  in  Figure  6d  turned  southward,  and  the  nose  of  the  magnetosphere,  as  determined 
by  equating  the  solar  wind  pressure  \n,yomj,v\^  to  the  magnetic  field  pressure  at 

the  subsolar  point,  moved  inweird  to  about  6.5  Re^  Here,  and  are  the  number 
density  (assumed  to  consist  of  H"'’  alone)  and  the  speed  of  the  solar  wind,  respectively, 
taken  from  Figure  7;  and  Bd  is  the  dipole  magnetic  field  strength  at  the  subsolar  point, 
taking  the  dipole  tilt  angle  at  that  time  into  account.  The  nose  position,  based  on  a 
smoothing  of  the  available  solar- wind  data,  for  25-26  August  is  shown  in  Figure  9. 
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Figure  6a.  IMP-8  IMF  data.  (Courtesy  of  R.  Lepping) 

15 


Magnetic  Strength(nT) 


Magnetic  Strength(nT) 


Magnetic  Strength(nT) 


Figure  7.  IMP-8  Solar-Plasma  data  (Courtesy  of  A.  Lazarus) 
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Figure  8a.  CRRES  magnetic  field  measurements  on  orbit  76.  (Courtesy  of  H.  Singer) 
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Figure  8b.  CRRES  magnetic  field  measurements  on  orbit  77.  (Courtesy  of  H.  Singer) 


MAGNETOPAUSE  POSITION  AT  SUBSOLAR  POINT:  AUG  25-27 


UNIVERSAL  TIME  (HR) 


Sharply  enhanced  fluxes  of  H'*',  He"^,  and  0"^  were  detected  with  IMS-LO  at  08:00 
UT,  when  CRRES  was  at  6>3  and  07:30  LTy  and  these  fluxes  were  repeatedly  observed 
in  the  region  3.2<L<6.8  and  03:00<MLT<10:30.  The  time  sequence  of  the  pitch-angle 
and  energy  distributions  of  these  ions  are  shown  in  Figures  3  and  4. 

3.S.2.2  Ion  Pressures 


In  order  to  describe  the  general  characteristics  of  the  ions  that  were  injected  and 
accelerated  in  the  trapping  region  of  the  magnetosphere,  the  parallel  and  perpendicdar 
pressures  of  the  ions  measured  with  IMS-LO  were  computed  for  the  duration  of  the  storm. 
When  the  IMS-HI  ion  fluxes  will  be  available,  their  pressures  too  will  be  computed  and 
added  to  the  low  energy  ion  pressures.  Since  most  of  the  partial  pressure  is  due  to  ions  in 
the  range  50  to  300  keV,  the  inclusion  of  the  higher  energ}*  pressures  will  provide  additional 
information  on  the  dynamics  of  the  magnetosphere  in  general  and  on  the  affect  of  the  ions 
on  the  geomagnetic  field  in  particular. 

The  partial  pressure  component  perpendicular  to  the  magnetic  field  is  given  by  the 
perpendicular  energy  component,  jmv^sin^a,  weighted  by  the  distribution  function,  /(v), 
and  integrated  over  velocity  space.  That  is. 


Pji-J  ^mv^sin}af{v)dv 

Sinularly,  the  parallel  pressure  is  given  by  the  equation 


In  terms  of  the  flux,  y{tc,a,T),  that  is  plotted  in  Figures  3  and  4,  the  pressures  as  a 
function  of  time,  T,  are  pven  by  the  equations 


Pj,{T)  =  27r 


;(u;,a,r)4tn*ada 


P|,(T)  =  4x 


f®ma« 

y/wdw  I  j{‘w,ayT)cos^asinada 


These  equations  were  integrated  from  Wmin  =  ^.11  keV  to  Wjnat  =  35  keV.  As  shown  in 
the  pitch-angle  distributions  of  the  ions  (Figures  3  emd  4),  the  measurements  are  nearly 
over  the  full  pitch-angle  range  a^in  >  10°  to  <  170°. 
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In  Figures  lOa-lOh  the  ion  pressures,  in  units  of  keV/cm*,  are  plotted  as  a  function 
of  time  (UT)  in  hours  for  orbits  75  though  82.  The  panels  starting  from  the  top  of  each 
figure  give  pressures  of  (1)  (2)  He"*",  (3)  O'^,  and  (4)  the  sum  of  H'*’,  He"^,  and  O'*". 

The  next  panel,  (5)  gives  the  plasma 


m  =  2mo 


PijT) 


which  is  the  ratio  of  the  total  perpendicular  pressure  of  the  ions  to  the  magnetic-field  energy 
density.  The  last  panel  (6),  gives  the  L-shell  of  the  measurements,  as  determined  from  the 
Olson-Pfitzer  magnetic  field  model.  In  the  upper  four  panels,  the  solid  and  broken-line 
curves  give  the  perpendicular  and  parallel  pressure  components  respectively. 

Note  that  during  the  orbit  preceding  the  storm,  t.e.  Orbit  75  of  Figure  10a,  nearly 
the  entire  pressure  is  due  to  and  that  its  perpendicular  pressure  exceeds  the  parallel 
pressure  at  all  times.  The  He'^  and  0*^  pressures  are  significantly  above  background  only 
at  low  L  values;  the  rapid  fluctuations  of  the  pressures  reflect  the  poor  statistics  of  the 
measurements.  Moreover,  the  plasma  0  is  small,  with  its  maximum  value  being  about 
0.6%  over  a  broad  region  at  L>5. 

On  the  next  orbit  (Figure  10b),  all  the  ion  pressures  increase  significantly  near  08:00 
UT.  Note  that  the  injected  is  very  nearly  isotropic  (the  pressure  components  are 
equal),  but  the  parallel  pressure  of  O’**  exceeds  the  perpendicular  pressure.  IVansport  cal¬ 
culations  reveal  that,  in  the  outer  magnetosphere,  fluxes  of  ions  that  recently  emerge  from 
the  ionosphere  are  anisotropic,  with  the  parallel  component  exceeding  the  perpendicular 
component.  The  maximum  value  of  0  of  the  low  energy  ions  is  about  3%  at  L:»6.1  during 
the  inbound  pass  of  the  satellite. 

On  orbit  77  (Figure  10c),  quasiperiodic  variations  of  the  particle  pressures,  with  pe¬ 
riods  of  about  one  hour,  occur  between  12:00  and  16:00  UT.  A  tentative  explanation  of 
these  variations  is  given  in  the  foUowing  section.  Note  also  in  this  figure  that  H'*'  is  still 
isotropic  and  that  the  He*^  and  0*^  have  parallel  pressure  components  that  often  exceed 
their  perpendicular  pressures.  As  revealed  by  the  0*^  pitch-angle  distributions,  the  lower 
energy  ions  at  high  L  values  generally  have  the  highest  pitch-angle  isotropies  (butterfly 
distributions).  These  too  are  discussed  in  the  following  section.  The  maximum  value  of 
in  this  orbit  is  about  2.3%. 

During  orbits  78  to  82,  the  H'*'  distributions  evolve  fairly  rapidly  toward  normal 
distributions,  with  the  perpendicular  pressure  exceeding  the  paredlel  pressure,  especially 
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Figure  10b.  Perpendicular  and  parallel  pressures  of  H+,  He"^,  and  0+  during  orbit  76. 
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Figure  10a.  Perpendicular  and  parallel  pressures  of  H'*’,  He'*’,  and  O'*’  during  orbit  75 
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at  the  lower  L  values.  The  He*^  and  0*^  fluxes,  however,  tend  to  remain  fairly  isotropic  at 
the  higher  L  values,  and  the  0  value  generally  decreases  during  the  recovery  phase  of  the 
storm. 


S.S.2.3  Interpretation  of  Data 


A  tentative  interpretation  of  the  IMS«L0  data  bued  on  the  observations  described  in 
Section  2.1  is  given  here.  Improved  understanding  of  the  pertinent  processes  is  expected 
after  all  the  CRRES  data  have  been  reduced,  crosB*correlated  with  CRRES  and  other 
(DMSP  and  Akebono)  data  sets,  and  discussed  with  the  various  experiments  and  theorists. 


S.S,2.3.1  Ion  Sources 

Both  the  southern  turning  of  the  IMF  and  the  compression  of  the  magnetosphere 
at  06:15  UT  on  orbit  75  indicate  an  enhancement  of  magnetospheric  convection  and  the 
presence  of  high  fluxes  of  high*density,  low>ene^  ions  moving  upward  along  magnetic  field 
lines  connected  to  the  cusp/cleft  ionosphere.  These  ions  are  composed  principally  of  , 
He*^,  and  0*^  that  are  accelerated  outward  from  the  ionosphere,  and  and  He'^'^  from 
the  magnetosheath  that  are  backscattered  and  degraded  in  the  upper  atmosphere.  Those 
ions  with  low  velocities,  near  and  above  the  gravitational  escape  velocity,  are  efficiently 
accelerated  by  the  enhanced  convection  field  and  transported,  within  one  to  two  hours,  to 
the  near>earth  plasma  sheet  in  the  nightside  nmgnetosphere  [Delcourt  et  al.,1989;  Cladis, 
1988]. 

The  increased  fluxes  of  He'*’,  and  O'*  were  detected  by  CRRES  at  Lss6.3,  07:30 
LT,  and  08:00  UT  (l|  hours  after  the  southern  turning  of  the  IMF).  The  pitch-angle  and 
energy  distributions  of  the  H'*  are  sinoilar  to  those  of  H**  in  the  plasma  sheet,  but  the 
distributions  of  He**  and  0'^  are  typical  of  the  ones  computed  for  ions  originating  in  the 
cusp/cleft  ionosphere.  Moreover,  the  observed  ions  were  dispersionless  within  the  time 
resolution  of  the  flux  measurements,  so  it  appears  that  the  ions  were  rapidly  transported 
inward  shortly  after  the  He'*'  and  O'*  reached  the  plasma  sheet. 

The  pitch-angle  distributions  of  the  O'*  measured  at  high  L  values  indicate  that  the 
locsd  ionosphere,  at  both  ends  of  the  field  lines,  is  an  additional  source  of  low-energy  ions 
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(<1  keV).  Note  that  the  lowest  energy  ion  have  extreme  “butterfly”  pitch-angle  distribu¬ 
tions  and  that  the  flux  increases  monotonically  from  pitch-angles  of  90*’toward  0°  or  180°. 
Toward  higher  energies  the  distributions  have  maxima  that  progressively  move  toward  90°. 
Hence,  it  appears  that  low  energy  0*^  originate  in  the  local  ionosphere  and  undergo  energy 
diiftision  as  they  move  upward  along  the  field  lines.  He*^  and  are  alto  expected  from 
this  ionospheric  source,  but  these  ions  do  not  have  pronounced  “butterfly”  distributions  at 
high  L  values.  A  possible  explanation  is  that  both  the  pitch-angle  and  energy  distributions 
of  upflowing  ions  are  modified  by  ion-streaming  instabilities. 

9.S.2.3.2  Variations  of  Particle  Pressures 

At  noted  earlier  the  four  humps  in  the  pressures  of  He*^,  and  0*^  shown  in 
Figure  10c,  on  orbit  77,  are  correlated.  Specifically  the  pressures  increase  and  decrease 
simultaneously,  implying  that  the  tame  process  is  affecting  all  the  ion  species  in  a  similar 
way.  Analysis  of  the  data  indicates  that  three  of  the  pressure  humps  are  due  to  changes 
in  the  compression  of  the  magnetosphere  and  one  is  due  to  an  additional  plasma  injection 
fi'om  the  magnetotail. 

The  first  pressure  rite  begins  at  about  12:20  UT.  Near  this  time  the  solar-wind  pres¬ 
sure  increased,  as  indicated  by  the  plasma  data  in  Figure  7,  such  that  the  nose  of  the 
magnetosphere  moved  inward  to  about  6  Re.  CRRES  was  in  the  nigbtside  magneto¬ 
sphere,  MLT%3-4  hr.  Since  the  inward  motion  of  the  nose  position  causes  ions  in  the 
nightside  magnetosphere  to  drift  outward  and  the  pressure  gradient  was  directed  inward, 
the  pressure  at  the  satellite  increased. 

The  next  pressure  rise,  which  begins  at  about  13:20  UT  appears  to  be  due  to  a  new 
injectioi  of  ions  from  the  tail.  The  increased  fluxes  were  detected  in  the  early  morning, 
MLT%4 .6-5.2  hr,  at  L%5-6,  beginning  about  one  hour  after  the  increased  compression  of 
the  magnetosphere.  Such  a  delay  of  an  hour  is  consistent  with  the  expected  transport  time 
of  ions  from  the  day  side  cusp/cleft  ionosphere  to  the  near-earth  magnetosphere.  Moreover, 
inspection  of  the  pitch-angle  and  energy  distributions  of  the  species  during  this  time  period, 
reveals  that  the  fluxes  increased  monotonically  toward  lower  energies,  resembling  typical 
spectra  of  new'ly  injected  ions. 

The  pressure  humps  at  14.7  hr  UT  and  16.2  hr  UT  appear  to  be  due  to  decreases  in 
the  compression  of  the  magnetosphere.  Based  on  the  solar-wind  data  shown  in  Figure  7, 
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the  nose  of  the  magnetosphere  moved  outward  to  about  8.3Re  at  14.7  hr  UT  and  to  10.9Re 
at  16.35  hr  UT.  At  these  times  CRRES  was  on  the  dawn  side  of  the  magnetosphere,  at 
MLT%6-7  hr  and  Ls«6.2-6.8.  In  this  region  ions  drift  outward  during  a  magnetospheric 
decompression.  Hence,  the  pressure  increases  since  the  pressure  gradient  is  directed  toward 
lower  L  values. 

The  large  depression  of  the  ion  pressure  centered  at  about  15.3  hr  UT  is  believed  to 
be  due  to  ion  precipitation  caused  by  an  mduced  electric  field  along  the  magnetic  field. 
This  is  a  new  loss  mechanism,  which  has  been  identified  during  the  study,  and  is  discussed 
in  the  following  section. 


d.S.2.3.S  A  New  Lots  Mechanism 

Inspection  of  the  pitch*angle  distributions  of  the  0*^  measured  on  Orbit  77  from  13:10 
UT  to  16:16  UT  reveals  that  the  flux  of  ions  less  than  a  few  keV  increases  steadily  as  the 
pitch'Sngle  increases  from  near  0*  to  near  180^.  Perhaps  the  most  extreme  example,  the 
one  measured  from  15:14  UT  to  15:45  UT  in  the  L  interval  6.67>6.81,  is  shown  in  Figure  5. 
Such  "clipped-wing  butterfly"  distributions  are  not  readily  explained  by  known  processes. 

In  Figure  11  are  shown  two  magnetic  field  lines,  projected  on  the  Y>Z  plane  (geomag¬ 
netic  solar-magnetospheric  (GMS)  coordinates),  that  bound  the  field  lines  traversed  by 
CRRES  during  these  measurements.  Note  that  the  satellite,  depicted  by  the  solid  sphere, 
was  in  the  northern  hemisphere,  so  the  ions  were  principally  directed  toward  the  equator. 
Therefore,  an  explanation  might  be  that  the  ions  were  somehow  absorbed  and/or  scattered 
downward,  along  fleld  lines,  by  instabilities  or  wave-particle  interactions  near  the  equator. 
However,  the  time  variations  of  the  pitch-angle  distributions  show  that  as  the  flux  at  high 
pitch-angles  increases  the  flux  at  low  pitch-angles  decreases,  and  vua  versa.  In  Figure  5, 
for  example,  the  fluxes  at  a^lbO**  of  ions  of  energy  greater  than  about  4  keV  are  imusu- 
ally  high  while  the  corresponding  fluxes  at  a^SO"  are  unusually  low.  The  ions  behave  as 
though  an  electric  field,  modulated  in  time,  is  present  along  the  magnetic  field  line  and 
directed  from  north  to  south. 

In  Figure  12  the  two  field  lines  discussed  above  are  projected  in  the  X-Y  plane,  still 
in  GSM  coordinates,  the  field  lines  appear  as  they  do  in  this  plan:  because  the  dipole  tilt 
angle  was  equal  to  about  20.5°  at  the  times  of  these  measurements.  Now,  note  in  Figure 
8b  that  the  B,  (actually  Bh)  component  of  the  magnetic  field  measured  on  CRRES  during 
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FIELD  LINE  TRACE  (TSYGANENKO)  1987  LONG  MODEL 
TIME=08:50:40,  DATE=17-MAY-91 
THE  X.Y.Z  COORDINATE  SYSTEM  IS;  MAGNETOSPHERIC 
KP=  >  5  .  EP0CH(Y:M:D:H:M:S)=90:  8:26;  15:22:40 

FIELD  LINE  IN  (Y,Z)  PLANE.  DIPOLE  TILT=  20.5  DEG. 
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Figure  11.  Magnetic  field  lines  on  which  "Clipped-wing  butterfly"  distributions 
were  observed.  In  CMS  coordinates. 
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this  period,  from  about  13.2  hr  UT  to  about  16  hr  UT,  was  fluctuating  appreciably  but  also 
generally  increasing  with  time  over  the  modeled  magnetic  fleld  strength.  After  converting 
the  magnetic  field  to  GSM  coordinates  and  assuming  Bz  to  increase  linearly  from  13.5  hr 
to  15.7  hr  UT,  then  dBtIdi  ss  7.6  x  10~*  nT/s  over  this  time  period.  If  this  rats  is  taken  to 
be  uniform  over  a  volume  extended  along  the  Zqsm  and  enclosed  by  the  field  line  on 
which  the  solid  sphere  (CRRES)  is  located  in  Figure  12,  then  the  induced  electric  potential 
over  that  field  line,  neglecting  the  presence  of  the  earth,  would  be  =  \AdBt/di\^  where  A 
is  the  area  enclosed  by  the  field  line  projected  on  the  X-Y  plane.  The  area  A  may  be  taken 
to  be  elliptical  with  semiaxis  a  and  b  equal  to  1.1  and  3.5Re,  respectively.  Then  A  =  wa6, 
and  V*  is  about  3.7  keV.  The  potential  drop  in  the  ionosphere  between  the  ends  of  the  field 
line  may  be  neglected.  So  the  potential  drop  along  the  field  line  from  the  northern  to  the 
southern  hemisphere  would  be  of  the  crder  of  4kV.  On  the  average,  the  observations  are 
consistent  with  a  potential  difference  of  about  IkV. 

The  pitch>angle  distributions  of  the  and  He*^  do  not  exhibit  such  pronounced 
“clipped'Wing”  distributions.  In  this  respect  the  pitch>angle  distributions  are  similar  to 
those  of  ions  of  difierent  species  that  flow  upward  from  the  ionosphere.  So  here  too  the 
pitch«angle  distributions  of  all  the  ^‘field  aligned*’  ion  species  may  be  affected  by  ion  stream- 
ing  instabilities. 

Although  the  proposed  mechanism  is  promising,  based  on  the  analysis  of  this  single 
event,  additional  events  of  this  type  must  be  examined  to  verify  the  mechanism.  The  data 
imply  ion  precipitation  in  the  southern  hemisphere,  so  we  have  been  seeking  measurements 
to  test  this  feature.  Using  the  Akebono  ephemeris,  which  we  acquired  from  NSSDC,  we 
found  that  that  satellite  was  very  near  an  appropriate  location  and  time  to  detect  the 
precipitating  ions.  It  was  moving  toward  the  south  pole  on  Day  238  1990,  and  at  16.5  hr 
UT,  it  was  at  L=6.58,  MLT=6.77  hr,  and  altitude=1868  km.  We  have  written  to  Mukai, 
inquiring  whether  he  obtained  measurements  at  that  time,  but  he  has  not  yet  replied.  If 
this  explanation  of  the  ‘’clipped- wing”  distributions  is  correct,  such  distributions  should 
occur  near  dawn  or  dusk  whenever  the  magnetosphere  is  being  compressed  or  decompressed 
while  the  dipole  tilt  angle  is  appreciable.  During  the  first  compression  of  this  magnetic 
storm,  the  dipole  tilt  angle  was  about  zero  and  “clipped-wing”  distributions  were  not 
observed. 


4  Future  plans 
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Agency  tape  processing,  the  generation  of  the  high  resolution  data  base,  the  survey 
plot  library,  and  the  summary  data  base  will  be  continued  on  the  SUN  on  a  routine  basis. 

Cross  calibration  of  IMS-LO  with  other  plasma  instruments  will  continue.  The  un¬ 
expected  revival  of  LOMICS  has  made  it  possible  follow  our  original  intention  to  cross 
calibrate  with  that  low  energy  mass  spectrometer.  However  we  will  continue  to  cross  cali¬ 
brate  with  LEPA.  Since  LEPA  does  not  separate  ion  species  it  will  be  necessary  to  select 
data  from  periods  which  are  dominated  by  a  single  ion  species  as  well  as  meeting  other  nec¬ 
essary  criteria.  We  will  also  attempt  to  cross  calibrate  with  MlCS  since  there  is  sufHcient 
overlap  between  the  energy  ranges  of  IMS-LO  and  MICS. 

The  preliminary  version  of  the  empirical  model  of  ring  current  ion  composition  will 
be  constructed. 

The  distribution  functions  of  low-energy  ions  (<35  keV)  measured  during  quiet  times 
will  be  used  to  estimate  the  quiet-time  electric  field  in  the  magnetosphere,  using  Liouville’s 
theorem  to  test  electric  field  models  along  ion  drift  paths  between  successive  CRRES  passes 
through  the  trapping  region.  As  time  permits,  this  analysis  will  progress  to  more  disturbed 
periods.  Further  examples  of  **clipped-wing”  pitch-angle  distributions  will  be  sought  and 
the  validity  of  the  new  loss  mechanism  described  in  section  3.3.2.3.3  will  be  investigated. 
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Appendix  1 


The  Low  Energy  Ion  Mass  Spectrometer  (ONR  307-8-1,2) 


41 


The  Low  Energy  Ion  Mass  Spectrometer  (ONR-307-8>1,2)  on  CRRES 
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Scientific  Objectives 

The  Low  Energy  Ion  Mass  Spectrometer  (IMS-LO)  is  part  of  the  SPACERAD  [1]  payload 
on  the  CRRES  [2]  satellite.  IMS-LO  is  designed  to  measure  the  composition  of  plasmas 
which  are  the  sources  of  radiation  belt  particles,  and  to  provide  data  on  the  origin  and 
acceleration  processes  of  these  plasmas.  In  order  to  achieve  these  objectives,  the  instrument 
measures  energy  and  mass  spectra  covering  the  ranges  of  E/q  =  0.11  •  35  keV/e  and  M/q 
ftor.  \  •  32  AMU/e  with  good  coverage  of  pitch  angles  throughout  the  CRRES  orbit. 

Measurement  Techniques 

The  IMS*LO*l,2  instruments  on  CRRES  rely  upon  essentially  the  same  design  that 
was  successfully  implemented  with  previous  Lockheed  instruments  on  SCATHA  (launched 
January  1979)  and  S3*3  (launched  July  1976).  Improvements  have  been  made  with  each 
version  of  the  instrument  in  range  of  coverage,  resolution,  and  operating  flexibility. 

IMS>LO-l  and  IMS*LO>2  are  identical  instruments  which  are  mounted  on  the  aft  surface 
of  the  spacecraft  with  their  look  directions  at  135**  and  105**  to  the  spacecraft  -l-Z  axis. 
In  the  spacecraft  X*Y  plane  they  are  both  at  45**  from  the  +X  axis,  clockwise  viewing 
in  the  +Z  direction.  CRRES  spins  about  the  Z  axis  with  +Z  oriented  sunward.  The 
look  directions  of  the  instruments  were  chosen  \.q  to  maximise  coverage  of  fluxes  near 
the  magnetic  field  line  direction  thoughout  the  mission.  Each  instrument  performs  ion 
composition  measurements  in  the  energy  per  charge  range  0.11>35  keV/e  and  the  mass  per 
charge  range  1  •  32  AMU/e.  The  energy  range,  which  is  covered  by  45  energy  steps,  is 
broken  into  three  contiguous  parts  (j  s  0,2),  each  consisting  of  15  energy'  steps  (t  =  0, 14). 
The  mean  energy  per  ^arge  (Eij)  of  ions  on  the  step  of  the  part  is  given  by 

Eij(keV)  =  0.109  x  (1.14008)^‘'^“^’^  (1) 

The  three  parts  of  the  energ}’  coverage  are  sampled  in  parallel  by  three  separate  analyzer 
and  sensor  units  ("heads”).  At  the  completion  of  each  15  step  sequence,  the  background 
counting  rate  is  measured  for  each  sensor  head.  The  mass  range  (1-32  AMU/e)  is 
covered  by  32  steps.  Alternatively,  the  spectrometer  can  be  commanded  to  a  Heavy  Ion 
Mode  as  described  in  the  section  on  Operating  Modes  and  Control.  In  addition  to  ion 
measurements,  each  of  the  two  instruments  monitors  the  background  electron  flux  at  four 
fixed  energies.  The  electron  channels  are  described  in  a  separate  section  below. 

Ion  Optics 

Each  of  the  IMS-LO  mass  spectrometers  consists  of  three  analyzer  heads  which  measure 
ions  in  a  different  portion  of  the  E/q  range  from  0.11-35  keV/e.  One  of  the  analyzer  heads 
is  illustrated  in  Figure  1.  Each  analyzer  consists  of  four  sections,  a  collimator,  a  velocity 
filter,  an  energy  analyzer  and  a  channel  electron  multiplier  detector. 

Ions  enter  the  instrument  through  the  collimator  which  provides  an  acceptance  cone 
of  approximately  5®  full  width.  Following  coUimation,  the  ions  enter  the  crossed  electric 
and  magnetic  field  velocity  filter  (Wien  Filter),  which  acts  as  the  mass  analyzer  (MA).  Its 
magnetic  field  is  fixed  while  the  crossed  electric  field  is  varied  according  to  the  value  of  E/q 
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and  M/q  being  sampled.  The  fields  are  oriented  so  that  the  electric  and  magnetic  forces 
on  an  ion  are  in  opposition,  and  the  electric  field  is  chosen  so  that  there  is  no  net  force  on 
an  ion  of  the  desired  mass  and  energy,  allowing  it  to  pass  through  the  MA  in  a  straight 
trajectory.  The  MA  magnetic  fields  in  the  three  heads  are  nominally  445, 1211,  and  3304 
Gauss.  The  electric  field  plates  of  the  MA  are  driven  symmetrically  and  the  voltages  for  a 
given  mass  selection  scale  as  the  square  root  of  the  ion  energy  which  is  determined  by  the 
voltage  applied  to  the  electrostatic  analyzer.  The  maximum  MA  plate  potentials  in  the 
three  heads,  corresponding  to  the  minimum  ion  mass  and  the  maximum  energ}'  per  charge 
sampled  by  head,  are  nominally,  ±35,  ±253,  and  ±1850  V  for  the  normal  mass  range  and 
a  factor  of  four  lower  for  the  heavy  ion  mode.  All  three  heads  are  driven  by  a  common  high 
voltage  power  supply  and  a  divider  network  which  provides  the  different  ranges  required. 

Following  the  velocity  filter,  ions  are  bent  180”  in  a  hemispberic.-I  sector  electrostatic  an¬ 
alyzer  (EA)  which  is  angle  focusing  and  energy  dispersive.  As  with  the  MA,  the  EA  plates 
in  all  three  heads  are  driven  by  a  common  power  supply  and  a  divider  network.  The  ratio 
between  consecutive  values  of  the  logarithmically  spaced  energy  steps  is  approximately 
1.14,  and  there  are  15  energy  steps  assigned  to  each  head.  Thus  the  energy  per  charge 
being  simultaneously  sampled  by  the  three  heads  are  spaced  by  a  factor  of  approximately 

with  the  first  head  covering  from  0.11-0.68  keV/e,  the  second  from  0.78-4.9  keV/e,  and 
the  third  from  5.6-35  keV/e.  The  EA  plate  separation  in  the  high  energy  head  is  reduced 
from  that  of  the  other  two  heads  in  order  to  achieve  the  desired  energy  range  without  the 
need  for  excessively  high  voltages.  After  exiting  the  EA,  ions  are  post-accelerated  by  IkV, 
to  improve  detection  efficiency,  and  enter  the  channel  electron  multiplier  detector. 

The  stepping  of  the  EA  and  MA  is  coordinated  under  digital  control  in  a  number  of 
modes  as  described  in  the  following  section.  Each  sample  of  a  particular  mass-energy 
value  takes  64  ms.  There  is  a  12  ms  deadtime  at  the  beginning  of  each  sampling  period  to 
ensure  complete  settling  of  the  analyzer  power  supplies,  even  in  the  worst-case  transitions. 

In  addition  to  the  energy  steps  described  above,  there  is  an  additional  setting  in  which 
each  of  the  three  heads  measures  background  coimts.  The  background  measurement  is 
performed  by  disabling  the  MA  power  supply  and  setting  the  EA  power  sup^'y  to  step  4. 
As  the  MA  power  supply  decays  toward  zero,  the  analyzers  pass  continually  higher  masse.<;. 
Following  the  12  ms  deadtime,  the  background  channel  is  sampling  mass  per  charge  values 
well  above  those  of  ambient  ions  and  measures  the  instrument  background  response  only. 

Electron  Detectors 

Each  of  the  IMS-LO  instruments  contuns  a  set  of  four  broadband,  fixed-energy  electron 
detectors.  Figure  2,  which  are  used  to  monitor  electron  fluxes  between  50  eV  and  25  keV. 
Electrons  entering  one  of  these  instruments  are  collimated  to  a  5”  acceptance  cone  and  pass 
through  an  aperture  which  sets  the  geometric  factor.  The  electrons  are  then  bent  180”  in 
a  magnetic  analyzer  and  pass  through  an  exit  aperture  which  sets  ^£/E  before  entering 
the  channeltron  sensor.  The  ^E/E  for  each  detector  is  50%  and  the  central  energies  of  the 
electron  channels  on  IMS-LO-1  and  IMS-LO-2  are  interleaved,  providing  a  total  of  eight 
channels  with  central  energies  from  0.067  to  20  keV.  The  electron  detectors  accumulate 
for  512  ms,  providing  samples  approxim<  .ely  every  6®  of  spacecraft  spin. 
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Operating  Modes  and  Control 

There  are  two  basic  submodes  of  operation  of  the  mass  spectrometers;  SWEEP  and 
LOCK.  Each  of  these  submodes  takes  32.768  seconds  to  complete,  approximately  one 
spacecraft  rotation. 

The  LOCK  submode  maximizes  resolution  in  time  and  pitch-angle  by  concentrating  on 
selected  masses.  The  instrument  is  locked  on  the  selected  mass  per  charge  while  stepping 
in  energy.  Each  energ}*  step  and  background  are  sampled  for  64m6  and  the  three  sensor 
heads  together  acquire  a  45  point  energy  spectrum  in  1.024  seconds,  approximately  12**  of 
spacecraft  rotation.  In  the  LOCK  submode  32  consecutive  energy  spectra  are  acquired  at 
the  selected  mass  before  the  instrument  moves  into  the  next  submode. 

The  SWEEP  submode  examines  the  mass  spectrum  in  more  detail  but  with  reduced 
resolution  in  time  and  pitch-angle.  In  this  submode  the  instrument  acquires  a  mass  spec¬ 
trum  at  each  energy  step  by  sweeping  through  all  32  mass  steps  while  the  energy  is  held 
steady.  Each  mass  step  is  sampled  for  64ms  and  a  full  mass  spectrum  is  acquired  in  2.084 
seconds.  The  complete  mass-energy  spectrum,  32  mass  steps  at  each  of  45  energy  steps 
together  with  32  background  measurements,  is  acquired  in  32.768  seconds. 

Eight  submodes  are  combined  to  form  the  three  basic  modes  of  operation;  LOCK-ONLY, 
SWEEP-ONLY,  and  SWEEP-LOCK.  As  the  names  imply,  LOCK-ONLY  consists  of  LOCK 
submodes  only,  SWEEP-ONLY  is  composed  of  SWEEP  submodes,  and  SWEEP-LOCK 
alternates  between  SWEEP  and  LOCK.  In  either  the  LOCK-ONLY  or  SWEEP-LOCK 
modes,  the  mass  steps  used  in  the  LOCK  submode  are  taken  sequentially  from  four  values 
which  are  held  in  the  instrument  memory.  These  values  may  select  any  four  of  the  32 
M/q  steps  and  correspond  to  any  four  ion  species;  for  instance,  Be*^,  He'*’'^,  0*^;  or 
0*^,  H*^,  0*^.  IMS-LO-1  and  IMS-LO-2  are  commanded  independently  so  that  it  is 
possible  to  have  the  two  instruments  operating  in  different  modes  simultaneously.  The 
operating  mode  and  LOCK  mode  mass  steps  can  be  selected  by  conunand.  Initially,  both 
instruments  were  operated  in  the  SWEEP-SWEEP  mode  while  their  performance  was 
assessed  and  on  orbit  measurements  of  mass  peak  shapes  were  made.  Subsequently  both 
instruments  were  commanded  into  the  SWEEP-LOCK  mode  with  the  LOCK  mode  steps 
set  for  H'*’  and  O’*". 

In  addition  to  the  normal  range  of  mass  coverage,  the  instruir.ent  can  be  commanded 
into  a  Heavy  Ion  Mode  for  use  during  chemical  release  experiments.  In  this  mode,  the  MA 
and  EA  voltages  are  reduced  by  a  factor  of  six,  scaling  the  M/q  coverage  upward  and  the 
energy  downward  to  give  an  operating  range  of  5  -  250  AMU/e  and  0.02  -  5.8  keV.  During 
barium  releases  one  of  the  instruments  has  been  operated  in  the  Heavy  Ion  SWEEP-LOCK 
mode  to  search  for  Ba'*'.  When  lithium  releases  were  made  one  of  the  instruments  was 
operated  with  the  LOCK  mode  steps  set  to  for  Li'*'.  The  other  instrument  remained  in 
the  normal  mode  of  operation  to  continue  to  monitor  ambient  plasma. 

The  discriminator  threshold  for  each  detector  is  set  by  conunand  to  one  of  four  values. 
A  special  calibration  cycle  will  be  conunanded  periodically  which  locks  on  a  selected  mass 
step  and  sequentially  steps  through  the  thresholds  in  order  to  monitor  the  channel  electron 
multiplier  gain.  In  addition,  the  calibration  conunand  may  optionally  enable  pulses  at  a 
known  rate  to  monitor  performance  of  the  signal  processing  electronics.  The  channel 
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electron  multiplier  bias  is  command  selectable  to  one  of  four  values  to  compensate  for 
possible  channel  multiplier  gain  degradation  over  the  life  of  the  mission. 

Calibration 

A  good  understanding  of  the  characteristics  of  the  mass  spectrometers’  response  to 
ambient  ions  is  necessary  in  order  to  correctly  interpret  the  data  it  obtains,  to  optimize  its 
operating  modes  and  to  determine  the  correct  ratios  for  the  £A  and  MA  voltage  divider 
networks. 

Three  approaches  were  used  which  provided  complimentary  information  on  the  spec* 
trometers’  characteristics.  During  the  design  phase,  detailed  ray  tracing  studies  provided 
a  baseline  description  of  the  variation  of  the  mass  spectrometers’  response  to  ions’  energy, 
mass  and  approach  direction. 

The  flight  instruments  were  individually  calibrated  in  the  Lockheed  Palo  Alto  Research 
Laboratorj’  mass  spectrometer  calibration  facility  in  which  they  were  exposed  to  monoen* 
ergetic  beams  of  selected  ion  species.  The  detailed  energy-mass  response  of  each  spectrom¬ 
eter  head  was  investigated  at  a  number  of  representative  energies.  At  each  energy  an  ion 
beam  of  H'*’  or  Ne'*'  at  the  selected  energy  was  directed  at  the  instrument  and  the  £A 
and  MA  voltages  varied  in  small  increments  to  provide  a  dense  mesh  of  measurements  to 
characterize  the  response  function. 

Since  the  full  response  is  dependent  on  approach  angle  and  incidence  position,  as  well  as 
energy  and  mass,  it  was  not  practical  to  perform  the  huge  number  of  individual  measure¬ 
ments  that  would  be  needed  with  a  monoenergetic  beam  in  order  to  fully  characterize  all  of 
the  spectrometers’  45  energy  steps.  Instead,  during  an  on  orbit  calibration  procedure,  am¬ 
bient  ring  current  ions  were  utilized  as  a  multi-ion  isotropic  extended  energy  source  which 
could  be  relied  upon  to  All  the  fleld  of  view  of  the  spectrometer.  With  the  instruments 
in  SW££P  mode  mass  spectra  were  acqmred  at  all  energy  steps.  These  measurements 
provided  relative  response  functions  for  common  mass  species  at  all  energy  steps.  Figure  3 
shows  peak  shapes  obtained  in  this  way  near  the  middle  of  the  head’s  energy  range.  The 
dashed  lines  represent  inferred  peak  shapes  while  the  marks  show  on-orbit  data.  Mass 
resolution  in  this  type  of  instnunent  is  dependent  on  mass  and  energy.  The  ^M/M  is 
lowest  for  at  the  low  energy  of  each  head  and  increases  with  increasing  energy  and 
mass.  At  the  middle  of  the  energy  range  dM/M  =  0.15  for  H'*'  and  Alil/M  ~  0.35  for 
0*^.  In  the  Heavy  Ion  Mode  ^M/M  =  0.16  for  Li'*'  and  0.47  for  Ba*^. 

The  acceptance  angle  is  approximately  conical  with  a  full  width  at  50%  sensitivity  of 
5”.  The  30s  satellite  spin  period  results  in  a  movement  of  the  instrument  of  less  than  one 
degree  during  the  64ms  sample  time.  The  instrument  geometric  factors  are  4  x  10~*  cm^ 
ster  and  ^£/£  of  %  10%  for  the  two  lower  energy  heads.  The  narrower  £A  plate  separation 
in  the  high  energy  heads  results  in  a  geometry  factor  of  2  x  10~*  cm^  ster  and  ^£/£  of  % 
5%. 

The  characteristics  of  IMS-LO  are  similar  to  those  of  its  predecesors  on  S3-3  and 
SCATHA  but  it  differs  from  the  ion  mass  spectrometers  on  some  other  magnetospheric 
satellites  such  as  ISEE-1,  D£-l,  and  AMPT£/CC£.  The  sensitivity  of  these  Instruments 
[3]  is  substantially  constant  over  their  energy  range,  but  their  field  of  view  and  AE/E 
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increase  substantially  at  low  energies  whereas  the  sensitivity  of  IMS-LO  is  proportional  to 
energy  and  its  field  of  view  and  idE/E  are  constant. 

Data  Reduction 

The  IMS'LO  instruments  completed  their  initialisation  sequence  and  began  normal  op¬ 
erations  on  August  4^^  1990.  They  are  operating  well  and  returning  good  quality  data. 
Figure  4  displays  an  example  of  sunoimary  plot  data  from  one  of  the  spectrometers.  The 
summary  plots  provide  an  exellent  overview  of  a  whole  orbit  and  display  mass  and  energy 
spectra  and  pitch-angle  distributions  as  a  function  of  time.  For  these  plots  the  data  have 
been  averaged  over  262s  intervals.  The  upper  panel  shows  mass  spectra  from  the  highest 
energy  analyser.  This  contains  horisontal  bands  near  the  bottom  and  top  which  correspond 
to  and  O'**  respectively.  Two  faint  intermittent  bands  corresponding  to  He*^  and  He'^'*' 
can  also  be  distinguished.  The  next  two  panels  show  energy  spectra  of  H'*',  and  0''‘.  The 
lower  two  panels  display  the  pitch-angle  distributions  of  and  0*^  For  each  of  these  the 
flux  has  been  sorted  into  10*  pitch-angle  bins. 
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Figure.  2.  Schematic  of  a  single  magnetic  electron  analyser. 
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Low  Energy  Ring  Current  Ions  During  the  Storm  of  26  August 
1900:  Observations  by  IMS>LO  on  CRRES 


H  L  Collin.  J  M  Quinn,  and  J  B  Cladis  (Lockheed  Palo  Alto 
Research  Laboratory,  3251  Hanover  St.,  Palo  Alto,  CA  94304) 

During  the  storm  of  26  August  1990,  about  two  hours  after  a 
sudden  commencement  which  occurred  at  0547  UT,  enhanced 
fluxes  of  He"*"  and  O’*"  at  energies  between  110  eV  and  35 
keV  were  detected  by  the  Lockheed  Ion  Mass  Spectrometer  on  the 
CRRES  satellite  at  0803  MLT  and  Lss6.4.  These  enhanced  fluxes 
persisted  for  about  three  days.  The  orbit  of  CRRES  enabled 
measurements  of  the  ions  to  be  made  repeatedly  in  the  region  3.0 
<  L  <  6.5  and  0300  <  MLT  <  1030.  The  time  sequence  of  ion 
observations  shows  the  development  and  decay  of  the  low  energy 
ion  component  of  the  plasma  in  the  ring  current  region  and  show 
evidence  of  multiple  particle  injections  and  ion  acceleration. 
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Siorm-Time  Dynamics  Implied  by  CRRES 
Measurements  of  Ring>Current  Ions 
J.B.Cladis.  H.  L.  Collin,  and  W.  £.  FVancis,  (all  at  Lockheed 
Missiles  and  Space  Co.,  Palo  Alto,  California  94304) 

On  26  August,  1990,  a  sudden  commencement  occurred  at 
0543  UT  and  sharply  enhanced  fluxes  of  He*^,  and  0*^  ions 
were  detected  with  the  Lockheed  Ion  Mass  Spectrometers,  IMS- 
LO  and  1MS*HI,  on  the  CRRES  satellite  at  0800  UT,  Ls«  6.38, 
and  0803  MLT.  The  enhanced  fluxes  were  observed  for  several 
orbits  in  the  region  3.2  <  L  <  6.8  and  0300  <  MLT  <  1030. 
Multiple  particle  injections,  ejections,  and  field-aligned  accelera¬ 
tion  events  are  implied  by  the  data.  The  distribution  functions 
of  ions  that  contribute  importantly  to  the  energy  dennty  of  the 
ring  current  were  used  to  compute  the  ion  pressures  parallel  and 
perpendicular  to  the  magnetic  field,  as  well  as  the  plasma  0  (the 
ratio  of  the  perpendicular  pressure  to  the  magnetic-field  energy 
density).  These  ions  include  H*^,  He"**,  and  0"^  that  were  mea¬ 
sured  in  the  energy  intervals  0.11  keV  to  1.6  MeV,  0.11  to  400 
keV,  and  0.11  to  100  keV,  respectively,  over  nearly  the  entire 
pitch-angle  range.  Graphs  of  the  ion  pressures  against  L  for  4 
CRRES  orbits  will  be  compared  with  pre-storm  values.  These 
resxdts  also  will  be  correlated  with  available  geomagnetic  indices, 
on-board  measurements  of  the  magnetic  field,  and  the  position 
of  the  magnetopause  at  the  subsolar  point  based  on  the  measure¬ 
ments  of  the  solar-wind  pressure  by  the  IMP-8  spacecraft. 


54 


Real  tima  monitoring  of  the  geophysical  parameters  in  support  of  the  CRRES 
chemical  releases. 

S.  B.  Mende,  H.  h,  Collin,  W,  X.  Peterson, 

Lockheed  Palo  Alto  Research  Laboratories,  Palo  Alto  CA,  94304 

T.  A.  Frits, 

Los  Alamos  Alamos  National  Laboratories,  Los  Alamos  NM  67545 
G.  Haerendel, 

Max  Planck  Institue  fur  Extraterrestische  Physik,  Garching,  Germany  8046 
R.  R.  Anderson, 

0.  A.  Hardy,  K.  J.  Singer,  0.  S.  Air  Force  Geophysics  Laboratory,  Ranscom 
AFB,  MA  01731, 

J.  C.  Foster,  MZT  Haystack  Observatory,  Nestford  MA  01686 

0.  L.  Reaaoner,  NASA  George  C.  Marshall  Spaceflight  Center,  Huntsville,  A1 

35812. 


The  high  altitude  releases  of  the  CRRES  program  were  active  experiments 
whose  success  depended  strongly  on  the  ambient  magnetospherie  geophysical 
conditions.  During  the  January  and  February  campaign  the  magnetospherie 
conditions  were  monitored  in  real  time  prior  to  release  in  an  attempt  to 
predict  conditions  for  each  release.  The  small  barium  releases  were  designed 
to  test  the  interaction  of  the  ion  clouds  with  the  ambient  magnetospherie 
environment  during  conditions  of  widely  varying  ambient  cold  plasma 
densities.  For  these  releases  the  in  situ  cold  plasma  density  was  monitored 
through  the  real  time  observation  of  the  wave  receiver  on  the  CRRES 
satellite.  The  large  barium  clouds  were  intended  t  produce  aurorral 
enhancements  at  the  foot  of  the  magnetic  field  line,  in  addition  to  the 
knowledge  of  the  cold  plasma  density  at  the  satellite  these  experiments 
required  the  observation  of  the  auroras  at  the  foot  of  the  field  line, 
measurements  of  the  ionospheric  electric  fields  in  the  conjugate  ionosphere 
and  reading  of  the  magnetic  field  and  energetic  particle  data  at  the 
satellite.  The  auroral  conditions  were  observed  by  the  instrumented 
aircrafts  and  the  ionospheric  fields  were  measured  by  the  Millstone  Hill 
incoherent  scatter  radar.  During  these  releases  the  required  geophysical 
conditions  were  met  which  contributed  to  the  success  of  the  experiments.  The 
large  Li  releases  were  intended  to  precipitate  local  fluxes  of  energetic 
electrons  by  pitch  angle  diffusion  in  the  presence  of  the  artificially  the 
injected  the  cold  lithium  ions.  This  experiment  dictated  the  toughest 
conditions  because  it  required  low  cold  plasma  density,  high  electron  fluxes 
at  the  satellite  and  uniform  dim  auroras  at  the  conjugate  region.  During  the 
campaigns  excellent  sets  of  data  were  gathered  which  characterised  the 
aurora,  the  simultaneous  particle  and  field  measurement  at  the  satellite  and 
the  accompanying  convection  in  the  ionosphere  during  various  conditions. 
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Distribution  in  Mi^netotul  of  0*^  Ions 
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Distribution  in  Magnetotail  of  O*^  Ions  from 
Cusp/Cleft  Ionosphere:  A  Possible  Substorm  Trigger 


J.  B.  Cladis  and  W.  E.  Francis 

Lockheed  Palo  Alto  Research  Laboratory,  Palo  Alto,  CA  94304 


ABSTRACT 

The  transport  of  0*'’  ions  from  the  cusp/cleft  ionosphere  to  the  magnetotail  during 
highly  disturbed  times  was  determined  by  computing  the  guiding>center  trajectories  of  the 
ions  to  a  distance  of  6  Rf  from  the  ionosphere  and  the  full>motion  trajectories  at  later 
times.  Case  histories  were  tallied  in  3  planes  perpendicular  to  the  Xgsm  uds,  3  planes 
perpendicular  to  the  Vosm  axis,  and  in  the  center  plane  of  the  tail.  At  various  times 
relative  to  the  enhancement  of  the  convection  electric  field,  the  following  ion  properties 
were  constructed  from  the  case  histories:  number  density,  mean  energ>’,  energy  and  pitch* 
angle  distributions  of  the  flux,  and  the  ion  pressure  components  parallel  and  perpendicular 
to  the  magnetic  held.  It  was  found  that  after  about  1.7  hours  the  ion  hux  in  the  near* 
earth  magnetotail  increased  dramatically,  and  the  spectrum  hardened  —  much  as  observed 
during  periods  just  preceding  substorms.  This  increase  is  attributed  to  (i)  the  increase  in 
the  0*^  outhux  from  the  ionosphere,  (ii)  the  increased  energisation  of  the  ions  by  the 
convection  electric  held,  and  (iii)  ion  trapping,  which  generally  occurs  because  the  ion 
magnetic  moments  generally  increase  after  the  ions  hrst  cross  the  geotail  center  plane. 
Moreover,  the  parallel  pressure  of  the  ions  exceeds  the  energy  density  of  the  magnetic  held 
at  Xqsm  <  -8  Re.  Based  on  the  expected  alterations  of  the  magnetic  and  electric  helds 
in  response  to  this  0*^  pressure,  a  substorm  trigger  mechanism  is  suggested. 
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INTRODUCTION 


During  times  of  enhanced  magnetospheric  convection,  the  upward  flow  of  high-density 
low-velocity  0*^  ions  flrom  the  dayside  cusp/cleft  ionosphere  increases  significantly,  as  dis¬ 
cussed  by  Yau  et  al.,  [1986].  Those  ions  with  energies  greater  than  about  the  gravitational 
escape  energy  flow  along  tailward-convecting  polar  field  lines  and  are  highly  accelerated 
along  the  magnetic  field  by  the  convection  electric  field,  especially  in  regions  where  the 
curvature  of  their  guiding-center  trajectories  in  the  earth  frame  is  large  [Cladis,  1986).  In 
the  nightside  magnetosphere,  the  ions  cross  the  center  plane  of  the  magnetotail  at  various 
positions,  depending  principally  on  their  initial  positions  and  velocities  as  well  as  on  the 
magnetic  and  electric  field  configurations  [Cladis,  1988). 

Delcourt  et  al.  [1989]  computed  the  three-dimensional  distribution  of  these  ions  in  the 
lobes  and  nightside  magnetosphere  using  the  magnetic  field  model  of  Mead  and  Fairfield 
[1975]  and  the  initial  convection  pattern  of  Volland  [1978].  Their  results,  while  repre¬ 
sentative  of  typical  distributions  during  quiet  to  moderately  active  conditions,  seriously 
understate  the  0*^  flux  in  the  midnight  sector  during  disturbed  times.  The  difficulty  lies 
in  their  assumption  of  the  conservation  of  the  magnetic  moment  of  the  ions  in  the  mag¬ 
netotail.  Indeed,  this  assumption  may  be  justified  for  some  trajectories  computed  in  the 
convection  electric  field  model  of  Volland  [1978],  which  is  uniform  over  the  polar  cap. 
But  more  realistic  convection  models  reveal  that  that  the  electric  field  is  intensified  ap¬ 
preciably  in  just  the  region  where  the  cusp/cleft  ions  are  convected  [e.g.,  Heppner  and 
Maynard,  1987].  Therefore,  when  the  potential  across  the  polar  cap  is  high,  ions  with 
initial  energies  so  high  that  they  usually  move  far  tailward  along  lobe  field  lines  are  now* 
convected  more  rapidly  toward  the  central  plasma  sheet.  Hence,  they  cross  the  center 
plane,  where  the  field-line  curvature  is  large,  with  much  higher  energies.  Under  these  con¬ 
ditions  the  magnetic  moment  generally  is  not  conserved.  The  ions,  prior  to  crossing  the 
center  plane,  axt  highly  aligned  with  the  magnetic  field.  Hence,  the  “scattering”  associated 
with  the  magnetic-moment  breakdown  about  the  center  plane  generally  causes  the  ions  to 
be  trapped  as  shown  below  (in  Figure  2).  Ion  trapping  greatly  increases  the  flux,  not  only 
in  the  region  of  the  breakdown,  but  especially  in  regions  along  the  inward  drift  paths  of 
the  ions,  where  the  ion  confinement  space  decreases  and  velocity  space  expands  due  to  the 
continued  acceleration  of  the  ions  by  the  convection  electric  field. 
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Here,  we  compute  the  time-dependent  characteristics  of  the  0"^  ions  in  the  nightside 
magnetosphere  as  a  function  of  time  measured  from  the  commencement  of  an  enhanced 
convection  event.  We  use  the  Heppner- Maynard  [1987]  convection  pattern  for  {By)iMF  >  0 
in  the  northern  henusphere  and  the  Tsyganenko  (1987]  “short-tail”,  active-time  [Kj)  >  5“ ), 
magnetic  held  model,  which  is  considered  to  be  appropriate  to  30  Re.  The  guiding-center 
equations  of  motion  are  integrated  in  calculating  the  ion  trajectories  to  a  distance  of  about 
6  Re  from  the  earth.  Beyond  6  Re,  full-motion  trajectories  are  computed.  The  results 
indicate  that,  within  1.7  hours  following  an  enhanced-convection  event,  the  spectrum  of 
the  0'*'  hux  in  the  near-earth  magnetotail  becomes  quite  bard,  being  comparable  to  the 
flux  measured  just  prior  to  substorm  onsets  [Daglis,  1990;  Kistler  et  al.,  1990;  Mobius  et 
al.,1987].  Moreover,  the  pressure  of  the  O'*"  ions  in  the  region  -13  <  Xcms  <  -8  Re, 
0.5  <  Ygms  <  2  Re,  and  -2  <  Zqms  <  +2  Re,  becomes  very  large,  with  the  0"^  pressure 
component  parallel  to  the  magnetic  held  exceeding  the  energy  density  of  the  magnetic  held 
At  Xqms  <  -8  Re.  The  accelerated  ions  therefore  affect  the  dynamics  of  the  magnetotail. 
In  fact,  correlation  of  the  timing  and  location  of  the  0*^ pressure  in  the  magnetotail  with 
observed  features  of  substornos  indicate  that  the  0*^  ions  may  play  a  dominant  role  in  the 
triggering  of  the  substorm.  From  the  changes  of  the  magnetic  and  electric  fields  implied 
by  the  0*^  pressure,  a  possible  triggering  mechanism  is  suggested. 

TRANSPORT  CALCULATIONS 

The  model  of  the  0*^ source  is  based  on  measurements  of  the  upflowing  ions  made 
with  the  Energetic  Ion  Composition  Spectrometer  (BIOS)  and  the  Retarding  Ion  Mass 
Spectrometer  (RIMS)  on  the  DE-1  spacecraft  over  the  dayside  cusp/cleft  ionosphere.  The 
statistical  study  by  Yau  et  al.  [1986]  of  the  EICS  data  obtained  during  the  period  Septem¬ 
ber  1981  to  January  1983,  normalized  to  the  reference  altitude  of  1000  km,  reveals  that 
the  outflux  of  O'*’  in  the  energy  range  .01  to  1  keV  peaks  near  the  invariant  latitude  of 
78®  in  the  noon  sector  (09  to  15  hr  MLT),  with  the  principal  outflux  lying  betv/een  about 
72  to  80®.  From  their  graph  of  the  O'*'  outflux  versus  invariant  latitude  for  disturbed 
times,  we  estimate  the  mean  outflux  in  the  in\'ariant-latitude  interval  72  to  80®  to  be 
about  5  X  10®  The  energies  of  O’*"  from  this  region  of  the  polar  ionosphere  are 

low,  ranging  from  about  2  eV  to  >  100  eV,  with  the  transverse  distribution  characterized 
by  a  temperature  of  about  10  eV  [Shelley  et  el.,  1982;  Waite  et  al.  1985;  Moore  et  al, 
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1986].  Accordingly,  the  source  was  taken  to  be  Maxwellian  in  the  upper  hemisphere  with 
a  temperature  of  10  eV,  normalized  to  an  outward  flux,  F,,  along  the  magnetic  field  of 
5  X  10® for  energies  higher  than  Wi  =  7.5  eV  at  the  altitude  of  1000  km.  That  is, 

2  r*  r*/2 

F,  =  — r  /  wdw  I  f{w,a)2Trsinacosada  (1) 

^  Jv)i  Jo 

where  f{w,a)  is  the  distribution  function  and  m,u>,  and  a  are  the  0*^  mass,  energy,  and 
pitch  angle,  respectively.  The  ion  source  was  taken  to  be  constant  from  09  to  15  hr  MLT 
and  from  72  to  80°  invariant  latitude. 

The  statistical  study  of  Yau  et  al.  [1986]  also  revealed  that  the  0*^ outflow  rate  in¬ 
creases  exponentially  with  the  Kp  index,  with  the  rate  at  Kp  >  6  being  a  factor  of  30 
larger  than  at  Kp  s  0.  It  is  therefore  expected  that  the  0*^  outflux  will  be  at  its  statistical 
active-times  value  during  an  enhanced  convection  event. 

MAGNETOSPHERIC  MODELS 

The  magnetosphere  was  modeled  by  (1)  the  ^'short-tail”  version  of  the  magnetic  field 
model  of  Tsyganenko  [1987],  as  improved  and  coded  by  David  Stem  of  NASA/Goddard, 
irith  the  model  parameters  set  for  Kp  >  5"  and  tilt  angle  =  0°,  and  (2)  the  convection 
electric  field  model  of  Heppner  and  Maynard  [1987]  that  is  generally  appropriate  in  the 
northern  hemisphere  for  {By)jMF  >  0.  The  polar-cap  potential  drop,  A'i',  was  taken 
to  be  100  kV,  which  is  appropriate  for  Kp  %  5°  and  AE  »  600  according  to  Table  1  of 
Heppner  and  Maynard  [1987].  This  correlation  of  Kp  and  AE,  based  on  the  polar-cap 
potential  drop,  is  consistent  with  the  data  presented  by  Rostoker  [1991]  on  the  general 
correlation  of  these  indices.  According  to  Caan  et  al.  [1978]  the  AE  index,  averaged  over 
1777  substorm  periods,  had  a  maximum  value  of  about  440.  Therefore,  the  models  used 
here  are  appropriate  for  more  highly  disturbed  times.  Both  a  more  relaxed  magnetic-field 
configuration  and  a  lower  polar-cap  potential  drop  would  lead  to  lower  0*^  pressures  in 
tbe  magnetotail.  The  conditions  adopted  ^ere,  as  discussed  imder  Results,  are  similar  to 
the  conations  under  which  the  substorm  data  on  0*^  ions  were  taken  that  were  used  to 
test  the  computed  fluxes. 

The  equi-potential  contours  of  the  Heppner-Maynard  model  for  A9  =  76  kV  are 
shown  in  Figure  1  (Fig,  3  of  Heppner  and  Maynard  [1987]).  An  efficient  routine,  which 
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is  described  in  the  Appendix,  was  developed  to  map  the  electric  field  from  the  reference 
plane  in  the  ionosphere,  where  the  convection  potentials  are  specified,  to  the  ion  positions 
in  the  magnetosphere. 

An  example  of  the  ion  motion  in  these  magnetic  and  electric  field  models,  which  also 
demonstrates  the  importance  of  the  full-motion  calculation,  is  shown  in  Figure  2.  Panel  2a 
shows  the  guiding-center  motion  of  an  O’*"  ion  under  the  conditions  Kp=4",  4,  4”^  and 
=  76  kV,  projected  in  the  XasMiZasM  pl&ne  ).  The  initial  position  of  the  ion  was  at  TQ** 
MLAT,  15  hr  MLT,  and  altitude  of  1000  km,  and  its  initial  energy  and  pitch  angle  were 
10.4  eV  and  92.5”,  respectively.  As  shown  in  the  figure  the  ion  crosses  the  center  plane 
at  about  Xqsm  =-14.5  Re  and  directly  enters  the  earth’s  atmosphere  in  the  southern 
hemisphere.  In  Panel  2b  the  full  motion  (broken  line  curve)  of  this  ion  is  shown,  starting  at 
the  radial  distance  of  6  Re.  The  solid  line  is  the  locus  of  the  center  of  curvature  determined 
from  the  instantaneous  positions  and  velocities  of  the  ion.  Note  that  this  trajectory  is  the 
almost  identical  to  the  guiding-center  trajectory  until  the  ion  first  crosses  the  center  plane. 
Afterwards,  the  ion  becomes  trapped  as  discussed  in  the  introduction. 

ANALYSIS 

The  running  time  on  the  VAX  11/785  required  to  compute  the  full  motion  of  the  ions 
was  appreciable,  even  though  the  routine  used  to  map  the  elec*  ric  field  is  more  efficient  than 
the  ones  used  previously  (e.g.,Cladis  and  Francis,  1989;  Delcourt  et  al.,  1989).  Therefore, 
the  number  of  trajectories  was  highly  restricted.  The  source  region  was  divided  into  24 
sections,  each  1  hr  MLT  by  2”  MLAT  at  the  1000-km  altitude,  and  6  ions  trajectories  were 
computed  fiom  each  section.  Because  the  case  history  of  an  ion  is  affected  much  more  by 
its  initial  speed  than  its  pitch  angle,  the  initial  pitch  angle  of  each  ion  was  taken  to  be 
the  same,  viz.  105”.  The  rix  initial  speeds  were  taken  to  be  at  the  centers  of  constant 
portions,  AF«,  of  the  parallel  flux  given  by  the  distribution  function.  That  is,  each  initial 
speed,  V/  was  determined  by  integrating  Eq.  (1)  from  to  with 

=7.5  eV  and  =  oo,  and  setting  AF,  =  |(5  x  10*) 

The  ion  case  histories  were  tallied  in  six  planes  perpendicular  to  the  Xqsm 
Xqsm  —  *5,  -6.5,  -8,  -10,  -13,  and  -20  Re;  three  planes  perpendicular  to  the  Yqsm  at 
Ygsm  =  -I*  -5,  and  2  Re;  and  in  the  center  plane,  Zqms  =  0-  The  listing  for  each  ion,  j. 
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in  the  plane  included  the  index,  s,  identifying  the  source  element,  the  plane  c’-ossing  point 
time  (T,),  velocity  pitch  angle  (ay),  total  energy  (tuy),  magnetic 

field  (J?,y,Byy,£xy,jBy),  electric  field  (£y),  and  \Ej  x  By/B?|  =  v^y.  Contributions  from 
the  ion  source  in  the  southern  hemisphere  were  taJcen  into  account,  assuming  symmetry 
about  the  center  plane,  by  adding  an  identical  listing,  but  with  2y  =  -2y,  v,y  =  -v,y, 
and  B*y  =  —  JB*y,  for  the  conjugate  ion  j*.  The  bulk  properties  of  the  ions  were  averaged 
in  certain  areas  ^4*  that  were  selected  on  the  basis  of  the  distribution  of  the  crossing  points 
in  the  planes.  The  areas  Ak  were  about  equal  to  {BtfBk)A,^  where  B,  is  the  field  intensity 
at  the  center  of  the  source  area,  and  Bk  is  the  field  intensity  at  the  center  oi  Ak- 


The  mean  ion  fiux  normal  to  summed  over  both  directions,  at  the  time  Tn  is, 
approximately, 


(2) 


where  vxt;  is  the  velocity  component,  perpendicular  to  the  plane,  of  the  ion  that  originates 
from  the  source  element  i  and  crosses  At  at  T  <  Tn,  and  v»y  is  the  total  velocity  of  that  ion. 
This  equation  is  not  strictly  correct  because  the  number  of  ions  per  second,  viz.  AF,Ai, 
represented  by  a  single  ion  crossing  within  A*,  may  not  all  fall  within  Ak  if  A<  does  not 
map  entirely  into  Ai,-  This  situation  certainly  occurs  when  ions  cross  Au  near  its  boundary’. 
However,  this  effect  tends  to  be  offset  by  ions  that  cross  the  plane  at  corresponding  points 
at  the  other  side  of  that  boundaiy.  The  error  in  the  flux  -*  and  in  the  additional  ion 
properties  listed  below  -  is  negligible  if  the  area  Ak  is  located  within  a  region  in  which 
the  density  of  ion-crossing  points  is  uniform,  but  it  becomes  appreciable,  about  a  factor  of 
two,  if  Ak  is  located  in  a  region  where  the  density  of  the  ion-crossing  points  is  high  and 
drops-off  abruptly  to  zero  beyond  one  of  its  boundaries.  In  general,  the  gradients  of  the 
ion  properties,  as  depicted  by  averaged  values  in  adjacent  areas  (histograms),  are  too  high 
by  amounts  that  increase  with  the  magnitudes  of  the  gradients. 


The  additional  ion  properties  averaged  in  Ak  are  the  following: 


The  nixmber  density. 


2 

v,y 


6? 


The  energy, 


1  J 

v)k{Tn)  =  j 


where  the  siunmation  is  over  the  ions  from  all  the  source  elenaents  that  cross  Ak^XT  <Tn, 


The  energy  spectrum  of  the  onmidirectional  iivix, 


Awi 


ii 

A*- 


where  Nu  is  the  number  of  ions  from  source  element  i  that  cross  Ak  with  energies  within 
Lwi  at  T  <  7n> 


The  ion  prewure  parallel  to  the  magnetic  field, 


where  tn||  is  the  ion  energ}’  component  parallel  to  the  magnetic  field  (w||  s  ^m(v  •£/£)*). 
The  ion  pressure  perpendicular  to  the  magnetic  field 


Vi) 


where  wj.  is  the  ion  energy*  component  perpendicular  to  the  magnetic  field  in  the  £  x  B 
frame. 


And  the  parallel  and  perpendicular  ion  pressures  divided  by  the  energy’  density  of  the 
magnetic  field  (B’/2mo)i 


and 


fAk^  Vij  B?,- 


where  B  is  the  magnetic  field  intenuty  at  the  ion-crossing  point. 
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RESULTS 


The  crossing  points  of  the  0"^  ions  in  the  planes  Xgsm=-^^  Yqsm-  0.5,  and  Zgsm— 
0  are  shown  in  Figure  3.  As  discussed  previously,  the  distributions  are  symmetrical  about 
the  plane  Zqsm—  0.  In  panels  (a)  and  (b)  the  crossing  points  at  Zg$m>  5  are  due  to 
the  ions  in  the  lobes.  The  points  at  the  lower  values  of  Zqsm  are  due  to  trapped  ions 
that  repeatedly  cross  the  planes  as  they  boimce  along  held  lines  and  drift  toward  the  west. 
Panel  (c)  exhibits  the  most  unexpected  feature  of  the  distribution;  nearly  all  the  crossing 
points  in  the  center  plane,  Zqsm  -  0,  occur  at  Yqsm  >  0  and  the  0*^  ions  gain  so  much 
energy  that  they  all  drift  toward  the  west.  This  feature  is  due  to  both  the  displacement 
of  the  ions  along  the  convection  electric  field  during  the  acceleration  process  and  the  large 

•a  «« 

Yqsm  component  of  the  E  x  jB  drift,  as  shown  by  the  potential  contours  of  Figure  1.  Note 
in  panel  (a)  that  the  ions  in  the  lobes  cross  the  Xqsm  -  plane  at  increasing  values 
of  Yosm  toward  decreasing  values  of  Zqsm>  ^ur  previous  calculations,  as  well  as  those 
of  Delcourt  et  al.[1989]  and  others,  performed  with  convection  models  in  which  the  equi* 
potential  contours  in  the  polar  ionosphere  were  parallel  to  the  Xqsm  axis,  the  center>plane 
crossing  points  were  more  nearly  symmetrical  about  the  Yqsm  plane. 

Histogranos  of  the  number  density  of  the  0*^  ions  in  the  center  plane  at  the  times  1.7, 
2.0, 2.3,  and  2.6  hours  after  the  conomencement  of  an  enhanced  convection  event  are  shown 
in  Figure  4.  These  densities  were  averaged  over  the  AXqsm  intervals  of  the  histograms 
and  0  <  Yqsm^  +2.  At  T  =  1.7  hr  the  number  density  is  already  fairly  high,  being 
about  0.1  cmT^y  and  nearly  constant  at  Xgsm<  *0.5.  Toward  later  times  the  number 
density  increases  everywhere  in  the  tall,  but  especially  in  regions  toward  the  earth,  owing 
to  the  inward  drift  of  the  trapped  ions.  In  the  plane  JfG5M=^*10,  the  density  averaged  in 
^Zqsm  bins,  with  0  <  Yqsm^  -^-S,  decreases  to  half  the  center- plane  value  in  the  bin  2  < 
Zgsm'S:  3.  A  statistical  study  of  the  ion  properties  in  the  plasma  sheet  in  the  region  (-10 
>  Xgsm>  *23),  measured  with  ion  spectrometers  on  the  lSEE-1  spacecraft  in  the  range  .1 
to  16  heV/e,  was  performed  by  Lennartsson  and  Shelley  [1986].  According  to  this  study, 
the  average  0"^  density  in  this  region  is  about  .2  cm”*  during  highly  disturbed  times.  The 
densities  shown  in  Figure  4  are  in  agreement  with  this  value  in  the  region  (A'g5m<*10) 
where  they  can  be  compared. 

Histograms  of  the  O'*' energy,  displayed  at  the  same  times  and  averaged  in  the  same 


64 


bins  of  the  center  plane,  are  shown  in  Figure  5.  Since  the  faster  ions  arrive  earlier  in 
the  tail,  the  mean  ion  energies  are  high  at  T  =  1.7  hr.  The  distribution  does  not  change 
appreciably  at  later  times,  especially  at  Xqsm<  *10-  The  mean  O'*’  energies  reported 
by  Lennartsson  and  Shelley  [1986],  which  were  averaged  over  AXqsm  intervals  of  5  Rb, 
[TosmI  <10>  *nd  |^c5mI  <10,  were  approximately  constant,  near  5  keV,  from  Xqsm  — 
-23  all  the  way  to  Xqsm  <  *6*  So,  while  the  computed  energies  are  is  agreement  with  these 
measurements  at  Xqsm  <*13,  they  are  appreciably  higher  at  Xqsm  >  *13.  Surely  some 
of  this  discrepancy  can  be  attributed  to  the  high  localization  of  these  cusp/cleft  ions  in 
space  and  time,  as  well  as  to  the  limited  energy  response  (<  16  keV)  of  the  fpecftroroeter. 
However,  as  discussed  in  the  following  section,  the  computed  energy,  averaged  in  the 
interval  .1  to  16  keV,  is  expected  to  be  lower  because  of  the  changes  in  the  magnetic  and 
electric  fields  that  occur  when  the  O'**  pressure  is  applied.  The  only  other  measurements 
of  0*^  ions  in  the  tail  were  made  with  the  Suprathermal  Energetic  Ion  Charge  Analyser 
(SULEICA)  on  the  AM?  f  ’ '  'IRM  spacecraft  at  A’’osm>*19  and  the  Charge  Energy  Mass 
Spectrometer  (CHEM)  ri  the  AMPTE/CCE  spacecraft  ai  XasM>'^'  High  fluxes  of 
O'*’ ions  with  energies  e^;  reeding  50  keV  were  measured  by  both  spacecraft  just  before-  and 
after-  the  expansion  phase  of  substorms  [Kistler  et  al.,1990;  Daglis  et  al.,  1990;  Mobius  et 
al.,  1987). 

The  energy  spectrum  of  the  omnidirectional  flux  in  the  center  plane,  averaged  over 
-5>ArcsAf^*13  and  0<yG5M<+2,  at  T  =  3.0  hr  is  shown  in  Figure  6.  Measurements 
of  ion  distributions  in  the  central  plasma  sheet  generally  indicate  that  the  distributions 
are  approximately  isotropic,  and  approximately  Maxwellian  at  energies  lower  than  that 
at  the  ‘‘knee”,  which  is  the  energy  above  which  the  directional  flux  decreases  sharply  to¬ 
ward  a  nonthermal  power-law  tail  [see,  e.g.  Christon  et  al.,1991).  Although  this  computed 
O’** distribution  is  not  the  usual  one  present  in  the  central  plasma  sheet,  the  energy  spec¬ 
trum  of  a  Maxwellian  with  a  “temperature”  oS  7.3  keV  (broken-line  curve)  is  drawn  in 
Figure  6  for  comparison.  Surely,  the  calculated  flux  is  not  sufficiently  detailed  to  assess 
such  a  feature  of  the  spectrum.  However,  the  pitch-angle  distribution  of  the  ions  is  clearly 
anisotropic.  It  is  a  butterfly  distribution  with  maxima  at  30  and  150®  that  are  higher  than 
the  90®  level  by  more  than  a  factor  of  4.  This  distribution  is  similar  to  the  one  measured 
with  the  AMPTE/CCE  spacecraft  near  Xqsm  =  *6.9  about  23  minutes  prior  to  the  occur¬ 
rence  of  a  substorm  [Daglis  et  al.,  1990].  Moreover,  after  dividing  the  flux  shown  in  Figure 
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6  by  a  reasonable  value  of  the  solid  angle,  fs  10  steradians,  it  is  in  satisfactory  agreement, 
over  the  entire  energy  range  (1  to  50  keV),  with  the  prc-substonn  flux  at  energies  >  1 
keV  reported  by  Daglis  et  al.  [1990],  and  with  the  pre-  and  post-substonn  fluxes  reported 
by  Kistler  et  al.  [1990]  at  energies  above  35  keV.  The  Kp  and  AE  indices  at  the  times  of 
these  measurements  were  similar  to  those  used  for  the  selection  of  magnetospheric  models. 
The  three-hourly  Kp  index  at  the  time  of  the  measurements  of  Daglis  et  al.  [1990]  was 
5'''  compared  to  the  model  parameter  5°;  and  the  AE  index,  averaged  over  the  pre-  and 
post-substorm  AE  values  of  the  12  substorms  included  in  the  study  of  Kistler  et  al.  [1990], 
was  633  compared  to  the  model  parameter  600. 

In  Figure  7  histograms  of  the  parallel  and  perpendicular  ion  pressures  in  the  center 
plane,  in  units  of  JbeV/em’,  are  shown  at  the  same  time  and  locations  as  those  of  the 
density  and  energy  in  Figures  4  and  5.  Note  that  the  parallel  pressure  (solid  lines)  exceeds 
the  perpendicular  pressure  (broken  lines)  at  all  times,  especially  at  Xqsm  <  *13.  From 
1.7  to  2.0  hours,  the  pressures  increase  rapidly  at  Xqsm  <  ‘6>3;  u  time  progresses  the 
pressures  increase  more  rapidly  in  successively  inward  bins.  Figure  8  shows  these  pressures 
divided  by  the  energy  density  of  the  magnetic  field.  Note  that  at  T  >  2  hr  both  of  these 
pressure  ratios  exceed  1  at  Xgsm  <  '8-  Since,  as  discussed  below,  the  perpendicular 
current  implied  by  these  ion  pressures  is  very  high,  these  0*^  ions  may  indeed  be  the 
current  carriers  sought  by  Kaufinann  [1987]  to  explain  the  pre-substorm  reduction  of  the 
magnetic  field  intensity  near  geosynchronous  altitudes. 

DISCUSSION 

The  O'*' pressure  shown  in  Figure  8  is  so  high  that  it  must  affect  the  dynamics  of  the 
near-earth  magnetotail.  Do  these  ions  trigger  substorms?  Phenomenologically,  it  appears 
that  they  do; 

(1)  They  are  injected  into  the  central  plasma  sheet  in  the  region  where  a  near-earth  neutral 
line  is  expected  to  form  [e.g.,  Nishida,  1984]. 

(2)  They  arrive  in  the  magnetotail  at  an  appropriate  time  (ss  1.7  hours)  following  the  com¬ 
mencement  of  an  enhanced  magnetospheric-convection  event,  considering  that  substorms 
often  occur  after  a  delay  of  1  to  2  hours  after  the  southern  turning  of  the  interplemetary 
magnetic  field  [e.g.,  Baker,  1986].  And, 


66 


(3)  Their  flux  is  consistent  with  the  pre-substorm  0"^  flux  measured  in  the  near-earth 
magnetotail  [Daglis  et  al.  1990;  Kistler  et  al.,  1990]. 

How  do  these  ions  trigger  a  substorm?  Bsdcer  et  al.  [1982]  suggest  that  the  presence 
of  0*^  ions  in  the  magnetotail,  owing  to  their  large  gyroradius  p,  may  trigger  a  substorm 
by  increasing  the  growth  rate  of  the  ion  tearing  mode  instability.  Their  suggestion  is  based 
on  the  derivation  by  Schindler  [1974]  of  the  linear  growth  rate  of  the  instability,  which  is 
proportional  to  Here,  another  mechanism  is  suggested  based  on  the  changes  of  the 
electric  and  magnetic  fields  resulting  from  the  arrival  of  the  ions. 


Clearly,  the  magnetic  field  configuration  must  begin  to  change  much  before  the 
0*^ pressure  exceeds  the  energy  density  of  the  magnetic  field.  The  change  of  the  magnetic 
field  in  response  to  the  particle  pressure  can  be  estimate  by  computing  the  perpendicular 
current  density  ji  in  the  center  plane  implied  by  the  0*^  pressure.  The  current  density  is 
related  to  the  particle  pressure  by  the  equation  [Parker,  1957], 


VPx  +  (Pii-Px)£ 


(3) 


where  e‘i  is  a  unit  vector  along  the  magnetic  field, « is  a  unit  vector  in  the  Xqsm  direction, 
B  is  the  magnetic-field  intensity,  and  Re  is  the  radius  of  curvature  of  the  magnetic  field. 
As  shown  in  Figure  7  at  T  >  2  hr,  P||  exceeds  by  more  than  a  factor  of  two  at 
XasM<  ’S*  Moreover,  Re*  which  is  less  than  about  1  Re  at  XosM<  *8*  is  much  less  than 
the  scale  length  of  VP±.  Hence,  the  second  term  in  (3),  which  gives  a  current  density 
that  is  directed  from  dawn-to-dusk,  is  dominant.  Taking  P||  -Px  =  2.5  IbeV/cm’,  Re  = 
1  Re,  and  B  =  9  nT  at  XesM  =  *9  where  VPx  ^  0»  then  j±  %  70  x  10“*®Amp/cm*. 
(The  dawn-to-dusk  current  density  at  this  point  in  the  magnetic  field  model  used  for  these 
transport  calculations  is  about  50  x  10”^®  Amp j cm}  [Tsyganenko,  1987].  Such  a  current 
density,  therefore,  causes  an  appreciable  elongation  of  the  local  magnetic  field. 

As  the  local  magnetic  field  elongates  under  the  pressure  of  the  ions,  an  electric 
field,  directed  from  dusk-to-dawn  in  the  center  plane,  is  induced  by  the  change  of  the 
magnetic  field.  Hence,  by  opposing  the  convection  electric  field,  this  induced  field  re¬ 
duces  the  pressure,  as  expected,  from  the  values  described  above.  The  induced  field, 
which  is  highly  structured  along  the  Yqsm  coordinate,  also  produces  velocity  shears 
in  the  earthward-convecting  flow  near  the  dawn  and  dusk  “boundaries”  of  the  injected 
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O’*" population.  Short  of  a  full  computer  simulation  of  the  dynamical  interactions  of  the 
plasma  and  the  magnetic  and  electric  fields,  it  is  difficult  to  assess  the  magnitudes  of  these 
effects.  However,  an  estimate  of  the  maximum  change  of  the  magnetic  field  during  the 
pressure  build-up  time  indicates  that  these  shears  may  be  very  large.  A  substorm  may, 
thv..ifore,  be  triggered  by  velocity-shear  instabilities  that  produce  short-scale  electric  and 
magnetic  turbulence.  Such  turbulence,  by  impeding  the  dawn-to-dusk  current  flow  across 
the  injected  0"^  region,  may  divert  current  to  the  ionosphere,  forming  the  current  wedge 
(e.g.,  McPherron  et  al.,  1973].  According  to  the  substorm-simulation  results  of  Kauf- 
mann  [1987]  and  Robinson  and  Rinaldi  [1991],  the  current  diversion  alone  could  trigger 
a  substorm.  Additionally,  however,  the  turbulence  may  enhance  field-line  merging,  which 
ultimately  provides  the  substorm  energy  [Stem,  1990;  Kan,  1990].  Velocity-shear  insta¬ 
bilities  in  the  auroral  ionosphere  are  discussed  by  Kintner  [1976]  and  Kelley  and  Carlson 
[1977].  Also,  the  Kelvin- Helmholtz  instability  has  been  applied  in  the  magnetosheath  by 
Southwood  [1979]  and  in  the  magnetotail  by  Rostoker  [1987]. 

CONCLUSIONS 

The  results  listed  below  strongly  indicate  that,  during  times  uf  enhanced  convection, 
the  O'*’  ions,  which  originate  from  the  cusp/cleft  ionosphere,  may  trigger  substorms; 

1.  In  the  near-earth  plasma  sheet,  the  flux  of  O'*'  increases  dramatically,  and  the  spec¬ 
trum  becomes  hardened,  resembling  the  O'*'  fluxes  measured  with  the  AMPTE-IRM/CCE 
spacecraft  at  periods  proceeding  substorms.  This  increased  flux  at  disturbed  times  is  due 
to  (i)  the  increased  outfiux  of  the  ions  &om  the  ionosphere,  (ii)  the  increased  energization 
of  the  ions  by  the  convection  eh  :tric  field,  and  (iii)  the  trapping  of  the  ions  resulting  from 
breakdown  of  their  magnetic  moments  during  their  first  crossing  of  the  center  plane. 

2.  The  increased  O'**  flux  and  pressure  occur  after  a  delay  of  1.7  hrs  following  the  en¬ 
hancement  of  the  convection  field.  Such  a  time  delay  is  consistent  with  the  observed  delay 
(1  to  2  hrs)  of  the  general  occurrence  of  substorms  following  the  southern  turning  of  the 
interplanetary  magnetic  field. 

3.  The  ion  pressure  is  applied  in  the  region  of  the  magnetotail  where  the  near-earth  neutral 
line  is  expected  to  form. 
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Estimates  of  the  magnetic  and  electric  field  chimges  in  response  to  the  O’*"  pressure 
in  the  tail  indicates  that  shears  in  the  convection  velocity  will  develop  near  the  flanks 
(dawn  and  dusk)  of  the  0*^  injection  region.  It  is  suggested  that  the  substorm  may  be 
triggered  by  velocity-shear  instabilities  that  effectively  create  anomalous  resistivity  in  the 
near  midnight  region  of  the  neutral  sheet. 

APPENDIX 

A  brief  description  is  given  here  of  the  routine  developed  to  map  the  electric  field  from  the 
ionosphere  to  the  ion  positions  in  the  magnetosphere.  Piecewise  mathematical  expressions 
of  the  convection  electric  field  distribution  given  by  the  potentials  are  pre-determined  as 
a  function  of  polar  coordinates  in  the  reference  plane,  and  the  corotation  electric  field  is 
added  vectorially,  giving  the  total  electric  field  (magnitude  and  direction)  as  a  function 
of  position.  The  routine  computes  two  mutually  perpendicular  electric  field  components, 
ejEsp  and  esEjp,  at  the  ion  position  by  tracing  only  two  magnetic  field  lines  as  follows: 
(i)  The  magnetic  field  line  through  the  ion  position  is  traced  to  its  intersection  (point  1) 
of  the  reference  plane;  at  this  point  the  total  electric  field  E  =  Eeonveetton  +  ^eorotation  >6 
computed  from  the  mathematical  expressions,  (ii)  A  second  magnetic  field  line,  displaced 
from  the  particle  position  by  the  distance  is  traced  to  the  reference  plane  (point  2), 
and  the  displacement,  rjSrfiy  of  point  2  from  point  1  is  determined.  Here,  17  is  a  unit  vector 
directed  from  point  1  to  point  2,  and  the  unit  vectors  ei,  ej,  and  es  are  in  the  directions  of 
the  magnetic  field  intensity  the  direction  E  x  VH,  and  the  direction  ej  x  62 ,  respectively. 
The  desired  components  are  then  obtained  from  the  equations, 

and 

where,  from  the  conservation  of  magnetic  flux, 

hi  _  Bj,  6(2 
6(3  Bi  6rj2 

Bp  and  Bi  are  the  magnetic  field  intensities  at  the  ion  position  eind  at  point  1 ,  respectively. 
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Computational  efficiency  is  further  improved  by  projecting  the  electric  field  forward 
for  several  time  steps  using  the  equations, 


where  i  =1,2  denotes  the  electric  field  components,  and  n  denotes  the  current  particle 
position.  The  exponents  are  determined  from  the  values  of  Ei  and  B  at  the  two  previous 
points  at  which  the  Ei  components  are  computed  by  tracing  field  lines.  The  final  projected 
values  of  Ei  are  compared  with  the  computed  values.  If  the  error  in  either  component  is 
excessive,  the  particle  is  returned  to  the  position  of  the  previously  computed  values  of  Ei 
and  the  number  of  projected  values  is  reduced  by  fraction  depending  on  the  magnitude  of 
the  error. 
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X(Re) 

Figure  2.  Guiding-center  trajectory  (panel  a.)  and  outer-magnetosphei.'e  portion  of  full- 
motion  trajectory  (panel  b.)  of  O'*'  ion  projected  in  Xgsm^Zgsm  plane.  Init‘  coudi- 
tions:  energy  =  10.4  eV,  pitch  angle  =  92.5®;  position  =  79®  MLAT,  1500  MLT,  and 
1000  km  altitude. 
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Figure  4.  Number  density  of  O'*"  ions  in  interval  0  <  Yqsm^  2  of  center  plane  versus 
XasM  times  shown  at  tops  of  graphs. 
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Rguie  5.  Mean  energy  of  0+  ions  in  interval  0  <  YasM<  2  of  center  plane  versus  Xasu  »t 
times  shown  at  tops  of  graphs. 
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Figure  6.  Energy  spectrum  of  O'*'  omnidirectional  flux  in  center  plane  averaged  over 
—5  ^Xqsm^  “d  0  ^YasM^  2  at  T  =  3.0  hr.  The  broken>line  curve  is  the 
spectrum  of  a  Maxwellian  with  a  “temperature”  of  7.3  kcV. 
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Figure  7.  Parallel  (solid  lines)  and  perpendicular  (broken  lii  es)  pressure  components  of 

0+  ions  in  interval  0  <  YasM<  2  of  center  plane  versus  Xqsm  at  times  shown  at  tops 
of  graphs. 
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PRESSURE  RATIOS,  0,  |  AND  0j. 


T  «  2.0  HR 


Figure  8.  Ratios  of  0*^  parallel  (solid  lines)  and  perpendicular  (broken  lines)  pressures 
divided  by  energy  density  of  magnetic  field  in  interval  0  <  Ycsm<  2  of  center  plane 
versus  Xqsm  at  times  shown  at  tops  of  graphs. 
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